A Statistical Description of Parametric Instabilities with an Incoherent Pump 
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The effect on parametric instability growth of pump wave incoherence is treated by deriving a set 
of equations governing the space-time evolution of the ensemble-average coupled-mode amplitudes 
and intensities. Particular attention is paid to establishing the regions of validity of the statis- 
tical description. Thresholds, growth rates, and amplification rates are given for both spatially 
and temporally incoherent pump waves. Both absolutely and convectively unstable modes are 
considered. The statistical results are verified where appropriate by numerical integration of the 
coupled-mode equations with different models of pump incoherence. 
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CLi; I. INTRODUCTION 

O ' The requirements of laser fusion targets for high power lasers with good laser beam uniformity has driven a quest 
' ^ • for new techniques for smoothing the intensity variations on the target surface. Early attempts at beam smoothing[l] 
J>-»' were not well characterized but more systematic techniques [2-6] have demonstrated significant improvements in beam 
uniformity. All techniques involve introduction of phase nonuniformities which replace the normal beam pattern, 
typically containing substantial hotspots, with a smaller-length-scale speckle pattern. Further addition of bandwidth 
to the laser provides temporal smoothing of the speckle. The primary motivation for investing in these smoothing 
schemes is to reduce the initial nonuniformities that can seed fluid instabilities such as Rayleigh- Taylor, yet there 
is also palpable interest in reducing the strength of laser plasma instabilities such as stimulated Raman or Brillouin 
scattering or two plasmon decay. Supplying a theoretical framework for understanding these laser plasma interactions 
with smoothed laser beams is the task we undertake in this article. In a subsequent article[7], the results obtained 
here will be applied to particular instabilities in geometries of interest. 

With the usual assumptions about slow variation of parameters with respect to the frequencies uoj and wave vectors 
kj of the wave, the coherent and incoherent problems can be studied within the context of the coupled mode equations 
(11.10 or II. 7). The wave group velocities and damping rates are denoted by Vj and Vj respectively; the strength of the 
coupling between the waves, 70, proportional to the pump or laser wave amplitude, has the dimensions of frequency 
and is the rate at which the waves grow without damping in an infinite homogeneous plasma. As the reader will 
remember, for a coherent driver, both convective and absolute instability may occur provided the laser intensity, i.e. 
, 7q, exceeds certain threshold values set by losses [8]. In an unbounded plasma, the convective coherent threshold is 

7o = ^1^2- (/■!) 

Absolute instability requires that the decay waves be oppositely directed (V1V2 < 0), and, the absolute coherent 
threshold is 



Temporal and spatial incoherence introduce two additional parameters, the temporal bandwidth Ao^o and the spatial 
bandwidth (or spread in wavevectors) A/cq which is inversely related to the correlation length AT^. 

The effects of both temporal and stationary spatial pump incoherence on parametric instabilities in homogeneous 
and inhomogeneous plasma has been studied extensively both theoretically and experimentally over the past thirty 
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years [9-68]. The first approach was quite naturally to consider the purely temporal problem with coupled mode 
equations in homogeneous plasma. Zaslavskii and Zakharov [9] considered the generic undamped decay instability 
with methods developed earlier for studying the relaxation of two level molecular systems driven by an external field 
[10]. They found that the convective growth rate was reduced, in the limit 70 <C Au), from 70 to ^q/Alo. Both Valeo 
and Oberman [11] and Tamor [12] addressed this problem with different methods and obtained the same result. Tamor 
included a damping rate on the Langmuir wave coupled to an acoustic wave and found a dispersion relation where 
the damping rate for the Langmuir wave nup is replaced by Vp + 2Au). He noted reasonably that the bandwidth was 
unimportant unless it exceeded the damping rate I'p but made no comment on the fact that the bandwidth appeared 
asymmetrically in his dispersion relation; that is, the apparent damping of the acoustic wave was unaffected. In a 
series of articles [13-15], Thomson made a number of important contributions. First, he noted that for equal damping 
on both modes, the threshold was 70 = fAw if Aw u and incorrectly speculated that, for unequal damping, the 
threshold was 7g = Acji^ii^2/inin(i^i, 1^2)- This guess was based on the assumption that, for growth to occur in either 
mode, the average amplitude of both must grow. 

The correct solution was presented by Thomson[15] later using an exactly solvable model [69]. For the purely 
temporal problem, the ensemble-averaged mode amplitude equations showed these equations decouple and have 
different thresholds, i.e. 

7^ = Awo2^i, (/.3) 

and 

7^ = AW0V2 (^-4) 

Thomson noticed and explained the asymmetry as due to averaging over the rapidly varying phase of the mode 
amplitude. He obtained a more appropriate threshold by solving the equations for the ensemble averaged intensities. 
The solution showed that the lower of the thresholds (1-3,4) was appropriate. We show, as did Thomson, that the 
incoherent growth rate for the average intensity is 

line = 270/ Awo (7.5) 

provided y\ < line < v^. This answer is expected, given that the average amplitudes grow with rate 7Q/Aa;o. 
However, Thomson did not point out the interesting fact (found also by Laval et al.[16] in the space-time problem 
when jVi) = )V2l) that, if 7j„c > ^i^^^^ the incoherent growth rate is 

line = 47o/Aa;o, (/-S) 

that is twice the expected rate. We will comment more on this in Sees. V-VII. Thomson then considered the space 
time problem with the ensemble averaged mode amplitude equations and again found an asymmetry in the thresholds 
for the absolute instability when Awq >> Vj- Thomson naturally assumed that the appropriate intensity threshold, 
in analogy with the temporal problem, was the lower one. Later work showed this to be incorrect. 

Laval et al. reconsidered [70] the space -time problem for the evolution of the intensities both by using the Bourret 
approximation and, in the case Vq — 00, the exactly solvable Kubo- Anderson process (KAP). They showed that the 
intensity thresholds are in fact much lower than the amplitude thresholds when Auq » Vj and, in usual cases of 
interest (1^2 > ^21^1/^1), the incoherent absolute threshold equals the product of the bandwidth and the damping on 
the slow wave 

ll = Auoi^2- {1-7) 

In a later publication [17], Thomson applied the results to stimulated Raman forward scatter and incorrectly concluded, 
as we discuss subsequently, that small amounts of bandwidth, comparable to its growth rate, can suppress this slowly- 
growing instability. During this time period, the nonlinear evolution of instabilities driven by broadband pumps 
was investigated by using particle-in-cell (PIC) simulations [18-20]. These studies demonstrated a reduction in the 
stimulated Brillouin reflectivity from 60% to less than 10% with 5% bandwidth. Because these plasmas were strongly 
inhomogeneous, no reduction is expected for modest bandwidth. Aw ~ 70, as shown by Thomson [15]. However, 
with Auj/luq ~ .05, the line separations in these simulations are larger than an acoustic frequency and each line acts 
independently. Kruer et al. [19] also suggest that bandwidth in an inhomogeneous plasma will be ineffective unless 
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Auj > Vq/1 where I is the interaction length set by plasma inhomogeneity. This can be rephrased in terms of the 
correlation length, Xc = Vq/Alu <— I where I ~ 7o/'«'(bii^2|)^^^ and k! = dk£^{x) / dx for — ko — ki — ^2; k^ = 
is the condition for phase matching. Kruer's condition is difficult to satisfy for typical laser systems. We return 
to this subject in our discussion of spatial incoherence effects. The correct theory [15] for convective amplification 
in an inhomogeneous plasma showed that the growth rate was reduced but the ensemble-averaged mode amplitude 
convectivcly saturates at the same level as the coherently driven amplitude. (Thomson actually misstated his result, 
although his analysis was correct. Later work [45,49] that solved the coupled mode equations clarified the answer.) 

Estabrook et al. also considered, with PIC simulations, the effect of laser bandwidth on stimulated Raman 
scattering again in inhomogeneous plasma. Bandwidth was represented by a series of equally spaced laser lines. 
A reduction in the reflectivity was found when the line separation was greater the growth rate. Each line acted as an 
independent pump, and the intensity was low enough that the instability was not strongly saturated. Direct compar- 
isons between homogeneous plasma theory and PIC simulations of SRS in a plasma slab were made by Forslund et al. 
[21] with good agreement for the dependence of the SRS growth rate on the bandwidth of the frequency modulated 
laser. These authors pointed out that, Aw/70 > 2, despite a fourfold reduction in the growth rate, only a modest 
reduction in the power reflected was observed in the final state. A further increase of bandwidth to Alo/^q > 10 
brought the instability below threshold. Later PIC work [24] with multiple lines also showed good agreement with 
theory for the Raman backscatter growth rate when 70 / Aoj <C 1 at several combinations of density and laser intensity. 
These simulations also showed a monotonic decrease with bandwidth in the absorption into hot electrons due to SRS 
and in the SRS reflected power. Less effect on the forward SRS was observed, consistent with the theory we discuss 
subsequently. Other theoretical work on pump bandwidth effects concerned applications to specific targets [30] , and 
the application to induced spatial incoherence (ISI) [44]. The effects of both temporal and spatial bandwidth in a 
ID inhomogeneous plasma [45] for a phase mismatch, ko{x) — ki{x) — k2{x) = k^{x), varying linearly, k^^{x) = k'x, 
or quadratically, k/^{x) = k"x^/2, was treated by analytical and numerical methods. In reference [47], the Green's 
function formalism of Brissaud and Fritsch [69], used by other authors [15,16,38] was reformulated in the language of 
effective Hamiltonian matrices. In a second paper[48], this formalism was used to investigate the competition between 
temporal bandwidth and inhomogeneity for a variety of parametric processes, including two plasmon decay. Early 
experiments [25,26] that attempted to observe laser bandwidth effects utilized plasmas that were too nonuniform to 
expect observable effects. Later experiments used gas jet targets [27,29] or used microwave plasmas [28]. Clayton et 
al. and Giles et al. divided the laser power into two lines and found a striking reductions in reflected SBS power 
compared to a single line with the same total laser power. The lines were separated by much more than a growth 
rate. 

It was recognized early [31-33] that random spatial modulation of the phase mismatch, (j) = J dxkA{x), could 
also reduce parametric instability growth rates. Using methods similar to the temporal problem, Kaw et al. used 
the steady-state limit of the coupled-mode equations (11-10) to find a convective gain rate. 7Q/4A|t'iW2|, where 
A =< Sk"^ > Ic- Here Sk is the wavevector mismatch induced by the random variation in the plasma properties 
and Ic is the correlation length for the random process. Further work [34-38,40-42] on random fluctuations of the 
phase concerned its effect on the growth in inhomogeneous plasmas where, for example, the stabilization of absolute 
modes by linear gradients in the phase could be undone by these fluctuations. In this article, we concentrate on the 
effects in uniform plasmas. 

Beam smoothing techniques use not only temporally incoherent but also spatially incoherent pump waves. In the 
focal plane, the laser wave can be considered a sum of randomly phased plane waves with a spread in wavevectors, 
Afco. In the limit of no temporal bandwidth, the spatial interference pattern of the pump wave envelope is stationary 
and is related to the effect of stationary plasma turbulence on parametric instability treated by Williams et al.[39] 
Using a novel approach (unrelated to the methods used in the present analysis), they considered the threshold and 
growth rate for absolute modes in a finite system of length L. This analysis finds that for suflaciently large spatial 
bandwidth, Ako, the fastest growth rate is reduced from the coherent absolute rate, 

when Uj = to the spatially incoherent absolute rate 

7 = 7oln(| Afco I L)/ ] Afco I (1 I + I V2 ]), (7.9) 

(To obtain 1-9 from Eq. 52 of Williams et al., the reader is advised that the correct normalization for the growth rate 
is 7o-\/| V1V2 1/(1 14 I + I V2 I). Their Eq. 3 is incorrect.) The incoherent result in Eq. (1-9) must of course be less 
than the coherent one in Eq. (1-8). The unexpected feature of (1-9) is the length dependent logarithmic factor which 
is related to the fact that the rate (1-9) is not the average rate sought in our analysis but the largest rate expected 
in a system of length L. Except for the logarithmic factor, we recover the scaling of the result of Eq. (1-9). In the 
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present article, we derive average amplitude equations by using the Bourret approximation [70] and average intensity 
equations by using the so- called random phase approximation (RPA). [71-74] The major objective of this article is to 
provide a unified theory of spatial and temporal incoherence effects on parametric instabilities. The approximations 
and assumptions necessary to arrive at the RPA equations (V.32) and the Bourret equations (V.l 4) that form the 
basis for theoretical results in Sec. VI are carefully explored and systematically presented in Sec. Ill through V. A 
comprehensive set of results is presented for the threshold and growth rate of both absolute and convective instability 
in the incoherent limit. The domains of validity of the statistical approximations are explicitly noted and outside these 
domains the coherent results are shown to apply. Moreover the same analysis is done for the spatial amplification of 
convectivcly unstable waves. Finally, we present in Sec. VII numerical solutions of the fundamental set of equations 
(11-10) for particular models of incoherence that illustrate the meaning of the averaging procedures and verify the 
main results. 

The RPA dispersion relation (VI. 1) for an infinite system obtained from the RPA equations (V.32) form the basis 
for the analysis of growth rates and thresholds in Section VI. There are two parameters that play a role in the ensemble 
average equations that measure the effective bandwidth, Acoj,j = 1,2 where 

Acoj = Auo I (1 - VjJVo) \ + \Ako- Vj \ (7.10) 

In the general case of spatial and temporal incoherence, there can be a gap between the domain of validity for the 
RPA dispersion relation and the coherent one. The coherent domain is the whole region where either spectral width 
Auj is less than the corresponding damping rate or growth rate. In this intermediate region, the average amplitude 
dispersion relations (VI. 2) are valid, and the more unstable one agrees with the RPA dispersion relation within a 
factor of two. It is on this basis that we argue that the RPA dispersion relation can be used in the whole domain 
(denoted the incoherent domain) complementary to the coherent domain. 

In the remainder of Sec. VI, a description of parametric instabilities is given including the conditions for absolute 
and convective instability. Both early time behavior, which is dominated by convective instability, and long time 
behavior, which is dominated by absolute instability (if it exists) or by spatial amplification are considered. It is 
worth noting that, in the incoherent domain, the average amplitude dispersion relations (VI. 2) allow only convective 
solutions whereas the RPA dispersion relation (VI. 1) allows the possibility of absolute instability if a threshold can be 
exceeded. However, the simpler average amplitude equations do provide the convective mode threshold (VI. 5), growth 
rate (VI. 7), and spatial amplification factor (VI. 13) and (VI.16) (within better than a factor of two) in agreement 
with those obtained in Sec. VI using the RPA dispersion relation. 

Absolute instabilities are of particular practical interest because, above threshold, the only limit to their growth 
is the finite laser energy or other nonlinear effects. On the other hand, an absolute instability generally has a larger 
threshold to overcome losses to collisional or Landau damping than a convective instability. For a coherent laser wave, 
the threshold laser intensity is determined by setting the spatial growth rate 70/ | V1V2 [^^^ equal to one- half of the 
sum of the spatial loss rates jyj/ \ Vj \ for the decay modes (Eq. 1-2). In addition, as given by Eq. (1-8), the absolute 
instability growth rate is smaller than the growth rate 70 by approximately the ratio (| V2/V1 |)^/-^ in the usual case 
I V2/V1 \« 1. The exact formula is also given by Eq. (VI. 11). With incoherent laser beams, the threshold is still 
determined by setting the effective spatial growth rate equal to the sum of loss rates but now the spatial growth rate 
is Max(7o/2Ati;j | Vj |), that is, the maximum spatial growth rate for the average amplitude equations. The general 
formula is given by (VI. 9). Note in the special but interesting case that Awi = Aui2 = Aluq, temporal incoherence 

increases the absolute threshold if Ra = {"/q/ Aluq){^\ V1V2 |/Min | |) < 1 which exceeds the naive criterion 
by the square root of the group velocity ratio. We observe that the statement Ra < 1 is equivalent to requiring 
that the coherent spatial growth length | V1V2 \^^^ /70 exceed the larger coherence length | Vi/Aloi \. A rigorous 
application[46] of the Bers and Briggs criteria[75] arrives at this same condition as a necessary condition (VI. 10) for 
the incoherent limit to apply. Above the incoherent threshold, the incoherent absolute growth rate is roughly 70 / Awi , 
(the exact formula is given by Eq. (VI. 8)) i.e. approximately the same expression as the incoherent convective growth 
rate but with a different domain of validity. 

There is a limit where the effective bandwidths Auij are determined only by spatial incoherence. Then | Aujj/Vj \ = 
Ako is independent of the group velocity and the growth rate, 

7f„^l)=47oV|Afco|(|V^i| + |^2|), (/.II) 

agrees with the scaling found previously and given in Eq. (1-9). Thus our general threshold and growth rate formulae 

indeed recover the correct limits of purely spatial and purely temporal incoherence derived previously. The coherent 
convective threshold, always lower than that for absolute instability, is given by the requirement that the growth 
rate 70 be greater than the mean loss rate (i^ii'2)^''^- The same criterion for an incoherent pump applies if the laser 
intensity is reduced by the ratio of line widths factor Rc = v\l Aujx + V2I Auj2- The effective bandwidth Aujj must 
be larger than the corresponding damping rate Vj for each mode to increase the threshold. Otherwise the coherent 
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threshold apphes (I-l). Far above the incoherent threshold, the incoherent convective growth rate is reduced by the 
factor 7o/Min(Awj). Note that when Awj ~ Awg the convective and absolute incoherent growth rates are equal in 
the overlapping domain of validity. More exact expressions are given by (VI. 5) - (VI. 7). 

For many cases of practical interest Aw^ ~ Auq because the spread in laser wavenumbers is sufficiently small 
and/or the group velocity is small enough that | Ako • Vj \« Awq and either | Vj/Vo \« 1 or Vj/Vo < 0. Then 
temporal incoherence is the dominant stabilizing influence. However two special cases deserve mention. Examination 
of the expression (VI. 8b) for the absolute instability growth rate above threshold shows that if the maximum value 
of the coherence length Vj/AiVj occurs for the minimum group velocity, the growth rate is strongly reduced 7"''* = 
V270/V1ACJ2 in the usual case | V2 \<<\ Vi \. This case can occur if there is a large spread in wavcnumbcr accompanied 
by a weak temporal incoherence so that 7"*** = 7o / I ■ ^1 | • The transient or convective growth rate for the same 
case (in the incoherent limit) is larger, -yconv = 7o/ | Afco • V2 |. The other unusual case occurs in the case of temporal 
incoherence for forward scatter when Vi ^ Vq- Then Aoji = Acjo(l ~ 'l^i/'K)) << Ai^o so that the coherent convective 
growth still occurs even if Awq ~ 70 • The case of practical interest is stimulated Raman forward scatter. [22-24] 

The fate of absolutely unstable modes necessarily requires consideration of nonlinear effects but, in a bounded 
plasma, the maximum amplitude of convectively unstable waves may be found by computing the spatial amplification 
rate. In general there are two roots that are most easily obtained by using the average amplitude equations (VI. 2) 
but are also accessible from the average intensity equations (VI. 1). A bandwidth Awq > 70, i^j is usually sufficient to 
reduce the spatial growth rate. For moderate intensities, i.e., below any absolute threshold, the spatial gain coefficient 
K ^ Jq/AluqVi (where | Vi |>| V2 |) is reduced by (7o/Awo)(V2/yi)^/^ in the weak damping case. At higher intensities, 
a second root with gain coefficient k ^ ^fQ/AujQV2 occurs provided Acjq > 7o(| Vi/V2)^^^ and 7q > Awoi^2, that is, 
above the absolute mode threshold when V1V2 < 0. Therefore both a convective and absolute mode exist; and, 
interestingly, the convective mode has a spatial gain rate that exceeds that of the absolute mode. A more detailed 
discussion and exact formulae are given in Eqs. (VI. 12) - (VI. 18). 

In the general case, it is difficult to summarize the effects of incoherence on convective and absolute stabilities. 
However, it is useful to consider a common case of practical interest where | V1/V2 |>> 1 and Awi ~ Auj2 = Acjq- 
Then diagrams showing stability regions can be constructed with axes Aloq/jo and 1^2/1^0, i-e. moving towards 
more incoherence in one direction and towards more damping in the other. Figure VI. 1 is such a diagram for 
convectively unstable modes. There are four regions: coherently unstable, incoherently unstable with reduced growth 
rate, coherently stable, and incoherently stabilized. In Fig. VI. 2, the diagram for absolutely unstable modes is drawn 
with four analogous regions. This figure graphically shows that, in only a small region of parameter space does 
incoherence reduce the growth rate but not completely stabilize. Figure VI. 3 shows the different regions for spatial 
amplification. Finally, in Fig. VI. 4, an overall diagram is shown for all regions for convective, absolute, and spatial 
amplification. 

To this point, we have presented results for infinite homogeneous plasma. For a large enough system, these results 
are a good guide to the behavior in finite systems. Nonetheless, real plasmas are finite and it is well known that 
a threshold length is necessary for absolute instability. Moreover in Sec. VII, numerical solutions of the coupled 
mode equations (11.10) are presented in a necessarily finite system. Thus partly as a guide to the numerical solutions, 
we find the threshold length for absolute instability in the coherent and incoherent limit as obtained by solving Eq. 
(II- 10) or Eq. (V.32) as appropriate. It is an interesting feature that the normal modes in the slab, sinusoidal in 
the coherent limit, are exponential in the incoherent limit. The threshold length given by Eq. (VI. 30) increases as 
expected with bandwidth above the coherent threshold length (VI. 29). 

Several features of the statistical description that provoked further analysis were the "unexpected" factor of two 
that appeared in the growth rate for the ensemble-average intensity, the question of the validity of the RPA description 
in the intermediate domain, and the applicability of these results to analysis of experiments using beam smoothing 
techniques. These aspects were examined in Sec. VII by integrating directly the coupled mode equations (11.10) 
or (VII. 1) with particular choices of random processes to represent the pump wave incoherence. For the purely 
temporal case, an analytic solution for the distribution function of mode amplitude (Eq. (VII. 10) is obtained which 
is remarkably broad if the damping rates are negligible. In fact, the width of the Gaussian distribution is equal to 
the mean. Thus, higher powers of the amplitude, e.g., the intensity grow at faster rates than the mean which gives 
rise to factor of two mentioned earlier. We compare this distribution to a numerically generated one in Fig. 6 for 
the same parameters. On the other hand, with sufficient damping, this factor of two does not appear: that is, the 
average intensity and amplitude grow at corresponding rates. In this case, the distribution is strikingly narrowed as 
also shown in Fig. 6. 

In the space-time problem this factor of two also occurs for absolute instability driven by an incoherent pump when 
the decay wave group velocities are equal in magnitude and opposite in sign. We have verified our supposition by 
numerical integration that once again the distribution is broad for | Vi |=| V2 | but is much narrower if | Vi \»\ V2 \ 
as shown in Fig. 7. 

The validity of the RPA description in the intermediate domain was examined numerically in Sec. VII. B, by 
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considering the model of a spatially incoherent pump driving a pair of decay waves satisfying the intermediate domain 
inequalities (VII. 13). Although we did show that the RPA equations appear valid, we also discovered that the 
distribution of mode intensities (Fig. 8) is unusual in that it consists of a slowly decreasing tail on a distribution with 
a peak at nearly zero growth. 

We also show in Sec. VII. C that the standard model of an ISI beam, Eq. (VII. 15), which has both phase and 
intensity variation, can be treated as an incoherent pump wave provided the temporal and spatial bandwidth are 
large enough. Thus with the appropriate identification of experimental parameters with Aloq and A/cq, the formulae 
in Sec. VI can be applied to experiments. [51-65] 

A few remarks are in order regarding the derivation of the statistical equations that form the basics for the results 
outlined above. 

The analysis begins in Section II with the completely nonlinear coupled-mode equations (II. 2) appropriate for the 

case when the pump and decay waves are weakly coupled and weakly damped. In the linear analysis of this article, 
the pump wave is unaffected by the decay waves and the characteristic growth rate of the parametric instability is 
simply related to the pump amplitude at its mean wavevector (II. 5). Normalization of the decay wave amplitudes 
to the average pump wave energy yields the linearized coupled-mode equation (II. 7) in Fourier space. From these 
equations, the envelope equations (II. 9) in Fourier space or (11.10) in real space are obtained by expanding the mode 
amplitudes about the value at the mean wavevector. These envelope equations are used in Sections III-V to obtain 
equations for the ensemble- averaged mode amplitudes and intensities. 

If the pump wave has a distribution of wavevectors, then a given pair of decay waves will be frequency matched to 
only a portion of the pump wave spectrum. Thus one is naturally led to consider the frequency mismatches (IH-l) for 
a general triplet of wavevectors or the mismatches (III. 2) at the mean value of the decay wavevector. It is assumed 
that there exist triplets for which there is exact matching so that the mismatch near the mean wavevectors is small. 
The maximum value of the mismatches at the mean decay wavevector determines that the interaction is coherent if 
the mismatches are both small in the sense defined by III. 6. 

Two approaches to generating equations governing the ensemble averaged behavior of the instability arc used in 
this article. The first employs the Bourret approximation [70] to obtain the dispersion relation (IV. 5) and (IV. 7) for 
the stability of the ensemble-averaged mode amplitude < Uj >. Each equation for < Uj > is simple in that it does not 
involve the other but it does involve integrals over the spectral density of the pump wave and the frequency mismatch. 
In the incoherent limit or Markov limit defined by (IV. 17), the dispersion relation takes the particularly simple and 
well- known form given by Eq. (IV. 13) which states that the coupling between the waves is reduced by the ratio of 
the characteristic growth rate to the maximum mismatch Aj max for i = 1,2. For the one dimensional case with a 
pump spectrum that is Lorentzian, e.g. a Kubo-Anderson Process, the dispersion relations (IV. 22) are the familiar 
ones derived previously[15-16] and the mismatch Aj rnax is equivalent to the spectral widths Awj defined by (IV. 20) 
that are, within at most a factor of two, equal to the pump wave bandwidth Aluq (except for the special case of 
forward scatter where the Aa>i << AllI2)- Note in one dimension, the problem of spatial and temporal incoherence 
are not independent in fact so our restriction to one dimension must yield the results[16] of Laval et al. However in 
our numerical models, this connection is broken for convenience and tractability. 

The average amplitude equations (IV.5 and IV. 7) are somewhat unsatisfying because < ai > and < a2 > evolve 
independently. For the case of temporal incoherence alone, this untidiness was remedied by obtaining equations for 
the ensemble averaged intensities which are coupled and symmetric. Nevertheless, the result obtained for convective 
instabilities with only temporal incoherence is just that obtained by consideration of the most unstable average 
amplitude. 

In Section V.B, the same procedure can be followed to obtain Bourret equations for the ensemble-averaged spectral 
densities (V.25) that are related to the mode intensities by (V.2). As with the amplitude equations (IV.5 and 
(IV. 7), these equations involve integrals over the pump spectral density and the frequency mismatch; but, in addition, 
the integral contains the spectral density of the other decay wave at the pump shifted wavenumber. Two further 
approximations can be made: first, the Markov assumption valid provided that the mismatch is larger than the 
maximum of the growth rate and damping rates, i.e. Mm{Aimax) > Max(7o,^'j) removes the spatial and temporal 
growth rate from the integrand; second, the assumption that the decay wave spectral densities are slowly varying 
functions allows these densities to be evaluated at the mean wavenumber at exact matching. The first approximation 
yields the Random Phase Approximation (RPA) equations (V.28); and the subsidiary second approximation yields 
the set (V.32) with which the instability analysis is done in Section VI. As they must do, the RPA equations (V.28) 
are shown to reduce to the correct intensity equations derived using the Bourret approximation when the pump wave 
group velocity is infinite. For the purely temporal incoherence modeled by the Kubo-Anderson process[69], this latter 
result is exact. 

Our conclusions are presented in Sec. VIII. There, we also give some examples of current interest and consider the 
special case when the temporal bandwidth is larger than one of the mode frequencies. 
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II. THE COUPLED MODE EQUATIONS 

A. The Coupled Mode Equations In Fourier Space. 

In this article we limit our discussion to the coupling of three wave-packets which are both weakly coupled and 
weakly damped. Wc assume that each Fourier component Al, of the wave packet i (with i = 0, 1, 2)) behaves in 
lowest order as a normal mode; i.e., it oscillates in time like exp(— iw^t), where is the real part of the frequency 
= wl — iul which characterizes the wave-packet i and which satisfies the dispersion relation Di{Cl\, k) = 0. The 
assumption of weakly coupled and weakly damped modes can be written quantitatively as 

\4\»4,^ {11-^) 

where v]. is the linear damping increment of the Fourier component A\, and where 7 denotes the inverse of the 

characteristic time for nonlinear evolution of the coupled modes. In this limit the equations describing the coupling 
of the three wave packets are of first order in time. In Fourier space they have the following general form; 

[dt + + Al = J dkoyko,ko<«-kj* (II-'^c) 

where A]^^ is proportional to the Fourier component of the electric field of wave i; wave refers to the pump wave 
and waves 1 and 2 to the decay waves. The coupling constants Vko,ki are derived from the usual field expansion of 
the fluid or the Vlasov equations. 

The parametric approximation consists in neglecting the RHS of Eq. (II. 2a) and in taking = Aj^^exp(— iWkpO 
in the RHS of Eqs. II. 2b and II2c. In this limit the latter equations describe correctly the parametric coupling of 
waves in the usual decay regime (the so-called modified decay instability and the modulational instability [76] cannot 
be described by these equations since they correspond to coupled mode equations in which the second order partial 
derivative in time must be retained). 

In this article we derive the conditions for a reduction of the parametric growth due to the pump wave incoherence. 
The incoherence of the pump wave may be either temporal, or spatial, or both. In the case of purely temporal 
incoherence, the wave-numbers fco of the different components of the pump wave all have the same direction; if Aloq is 
the spectral width in frequency space of the pump wave and is small compared to the pump frequency, Awq << ujq, 
the spread AKq of the corresponding wave- number modulus fco =| fco | is given by 

Afco ~ Ai^o/Vgo (J/.3) 

Here Vgo is the characteristic group velocity of the pump wave and Vgo denotes Vgo = | Vgo | • In the case of purely 
spatial incoherence, all the wave- numbers ko have the same modulus fco =| ko | and are spread within a small cone, 
the half angle of aperture of the latter being denoted by AOq. In the general case, the mean wave-number and 
frequency of the pump wave are denoted by Kq and loq = Uj^^ respectively, and the incoherence of the pump wave is 
characterized by the spectral width Afco = oJo/Vgo and by the half-angle spread A6(). We will assume an azimuthal 
symmetry around Kq and a cross section of the total domain of existence of ko is displayed in Fig. 1. 

An assumption made implicitly when using Eqs. (II. 2) to investigate the incoherence effects is that the spread in 
frequency At^o is much smaller than the frequencies o;^. of the various waves. Similarly, the spread in wave number 
I Ako |~ Auio/Vgo + KqAQq is assumed to be small compared to wave numbers | kj |. 

We denote by 0:2 the angle between the propagation directions of the pump wave and of the wave 1. For a given 
direction 02 of the wave 1, there is a unique couple (Ki, w^J, (^2,^X2) satisfying the two resonance conditions, 



Ko = Ki+K2. {II A) 
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Throughout this article we investigate the reduction of the parametric instabihty corresponding to a given direction a2 
and the capital letters Ki will be used exclusively for this wavenumber of wave i which corresponds to the resonance 
conditions (11.4), for i = 1, 2; as stated before, Kq denotes the mean wavenumber of the pump wave, and the symbol 
ko will denote a generic wavenumber of the pump wave; lastly ai is the angle between the mean pump wave-number 
Kq and the wave number K2 corresponding to the resonance conditions II. 4 written for an angle 02. 

We may now relate the characteristic growth rate of the parametric instability corresponding to the coherent case, 
7o, for a given angle of observation a2, to the mode coupling constants and the pump amplitude by 

7o' =1 ^Ko,K, I' (I Ao |2) (77.5) 

where ([ Aq p) denotes the average energy of the pump wave in a sense defined in the next section. 
Defining the dimensionless coupling constants and wave amplitudes by 

^^ko.ki = Vko,ki/Vko,Ki (77.6a) 

«L = ^L/ <l ^0 |'>'/' {IIM) 
The dimensionless form of equations is simply 

[dt + iwk, + t-kj =^0 j (ikoi^ko.kiako (ako-kj* (77.7a) 



[dt + iul^ + fkj ak2 = To / <^koVko,ko-k20ko ("ko-ks)* (77.76) 

where I'ko.ka is a slowly varying function of ko and k^j, with uko.Kc = 1- Finally, the latter equations are supplemented 
by the following equation for the pump wave amplitude 

< = <exp(-i<t) (77.7c) 



B. The Envelope Approximation for the Coupled Mode Equations in Real Space 

At this point we can make a connection between the coupled mode equations (II. 7) written in their general form 

in Fourier space and their so-called envelope equation form in real space. The envelope approximation for the 
coupled mode equations corresponds to the first order expansion of w^. in a power series of (kj — K^). By writing 
Wjj, ~ wj^. + (kj — Kj) • Vgj, and by setting 

(i^i = «ki-Kiexp(-zwKi*) (77.8a) 

or i = 1 and 2, and 

«ko = flko-Koexp {-i[co°Ko + (ko - Ko) • V^o]*) (HM) 
the coupled mode equations can be written as 

{dt + iK ■ Vgi + 1^1) ak; = 70 /" dk(,a°, (a')k;,_k; ^xp - « (k(, • V^o) t (77.9) 



{dt + ik2 • Vg2 + 2^2) Sk' = 70 / dk^a^, (a^)* _ , exp - i (k^, • Vgo) t 

A J ^ ' ^0 2 

where k^ = kj — Kj and we neglected the slow variation of i^Ki+k' and of fKo+k^.Ki+k'^ with k.^ and kg. By taking 
the inverse Fourier transforms of Eqs. II. 9, one obtains the mode-coupling equations in their envelope approximation 
limit, namely 



{dt + Vgi ■ 9x + 1^1) ai(x, t) = 'joS (x - Vgo*) 
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[dt + Vg2 ■ + V2) a2(x, t) = 7o5 (x - Vgo^) al (77.10) 
where the quantity S'(x — 'Vgot) is defined by 

5(x - Vgoi) = j dk(,ak^exp(zk(, • (x - V^oi)) (77.11) 

Equations 11.10 arc the envelope equations which have been investigated with the statistical methods described 
in the next section. These equations are generalizations to three dimensions of one dimensional equations used in 
previous treatments of this subject[16]. For the sake of simplicity the envelope approximation will henceforth refer to 
either the coupled mode equations written in their envelope form 11.9, or the first order expansion of lo].,. 

Concerning our notation, Vga(ka) will denote, in the following the group velocity corresponding to a generic 
wavenmnber ka, namely Vga(ka) = {dw^^/d^ia); the notation Vga will be reserved to the case of exact resonance, 

i.e., Vga = Vga(Ka). 

On the other hand, the quantity (x,t) will denote the Fourier transform of the field a^^{t) in the case where 
one considers the mode coupling equations in their general form (II. 7); the notations aa{x,t) and €l0^ will be used 
by contrast in the case where the fast time and space variations have been factorized ab initio as in Eq. (II. 8a). The 
connection between the two sets is given by 

= aK„+kj,exp(+ia;^^f) (77.12a) 

and 

aa{x, t) = da{x, t)exp z(Kq • x — uj^h.) + c.c. (77.126) 



III. STATISTICAL DESCRIPTION FOR COUPLED MODE EQUATIONS 

A. The Frequency Mismatches 

In the two following sections, we derive equations describing the evolution of average quantities such as the average 
amplitude (oq), or the average intensity (| aa P) of wave a. Here < (f) > denotes the statistical average of the physical 
quantity (f), and its fiuctuation is written as S<i> = (j)— <((>>. The meaning of a statistical average can in some sense be 
understood as a time average, and the use of a statistical framework is justified for the following reasons: a reduction 
of the parametric growth can be expected whenever the spectral width || Ako || of the pump wave is large enough 
(due to a temporal or spatial incoherence) for its correlation time - as seen by the decay waves, in a sense to be defined 
in this section - to be shorter than the other characteristic times, namely, the damping and the growth time, and 
7^^. In this case, the statistical description is justified whenever the correlation length in K-space of the pump wave 
is small compared to its spectral width (e.g., when the number of ISI echelons or phase elements is large). In this 
limit the pump wave electric field can indeed be regarded as a stochastic process with a short correlation time and the 
standard statistical techniques can be applied. In the case of a purely temporal incoherence, the statistical description 
is physically justified only if the pump wave bandwidth is caused by some stochastic process, e.g., the spread over 
several independent lines of a non-monochromatic laser. Similarly, in the case of a purely spatial incoherence, the 
natural spread in angle of the pump wave wavenumbers must follow from the sum of many independent beamlets. 
Such a statistical independence may result from a random phase shift given to the different beamlets by means of 
a transparent phase mask [3]; it may also be due to the scattering of the incident pump wave upon static random 
fiuctuations. Lastly such a spatial incoherence can be achieved in the so-called ISI technique[2] by a combination of 
delay increments At given to the beamlets by echelon structures and of the laser temporal incoherence in the case 
AwoAt > 1. In these techniques, the laser beam is broken up into a number of statistically independent "beamlets" 
that are brought to focus by an optic of f- number, f^. In the focal plane these beamlets overlap creating a spatially 
nonuniform pattern with coherence length across the beam Xc = tt / (kosinO f) where 9f is the half-angle of the optic, 
i.e. 9f ^ (2/#)-i. 

We will henceforth restrict ourselves to physical situations where the statistical independence between the pump 
wave Fourier components is satisfied; the pump wave electric field will thus be regarded as a stochastic variable with 
zero mean; its statistical properties will be assumed to be entirely determined by its spectral density denoted as n'j^^ , 
the latter being itself characterized by its spectral width Ako, that is to say by Afco =| Ako |= AuiofVgo (temporal 
incoherence) and A^o (spatial incoherence). 
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Before deriving the statistical equations for the time evolution of < ai > or <| p>, let us first consider on 
a physical basis the reduction of the parametric growth caused by the pump wave incoherence. Suppose that we 
consider the parametric coupling corresponding to a given angle a2 between the pump wave and the first decay wave 
direction; this angle defines the resonant wavenumbers Ki and K2 of the decay waves corresponding to the exact 
resonance condition II. 3 with the mean pump wave number Kq. If one considers now a generic wavenumber ko 7^ Kq 
of the pump wave, the resonance conditions II. 3 can no longer be satisfied for the same wave-numbers Ki and K2. 
One is thus led to define the resonance mismatches Ai and A2 as 

Ai(ko, k2) ^ - a;W_k. - 4? 

and 

A2(ko, ki) ^ - 4;) - 4t^^ (7/7.16) 

When discussing the effect of the pump incoherence in the case of a given geometry defining the resonant wave 
numbers Ki and A'2: it is sufficient to consider the rjuantities Ai(ko;K2) and A2(ko,Ki). In order to simplify our 
notation, we will denote in our paper by 9 the value of a quantity 9 in the case of the exact resonance ki = Ki and 
k2 = K2. Accordingly we define Ai(ko) and A2(ko) by 

Ai(ko) = Ai(ko,K2) (7/7.2a) 

A2(ko) = A2(ko,Ki) (777.26) 

In the case of an incoherent pump wave, the spectral width Ako gives rise to a characteristic size for Ai(ko), which 
will be denoted as Ajmax, i-e. 

Ai „,a. = A(Ko + Ako) {111-3) 

For instance, in the case where the envelope approximation can be applied, one readily finds Ai max = (Vgo — 
Vgi) • Ako which can also be written as 

Ai max =1 {Vgo - Vgicos a2)Ako - TCoV^gisin Q!2Aeo | {IIIAa) 

or 

Aimax =1 (1 - Vgicos a2/Vgo)Awo - -Fl'oV^isin a2Aeo | {III Ah) 

More generally, it has been shown[7] that it is sufficient to expand Ai max as follows: 

Ai„ax =1 ii Afco + ai AOo + /3iAeg | (777.5) 

with t\ = {Vgo — V^iC0Sq;2), fJi = —VgxKosm.a2- For backscatter, a\ vanishes and it is necessary to expand to 
second order in A0o. The coefficient of AQq,/3i, can be computed in terms of the waves parameters duj" / d'K.j , and 
9^cj"/9Ki9Kj, and the angles a-i defining the geometry of the interaction. The first term in (III. 5) corresponds to 
temporal incoherence and the last two to spatial incoherence. 

In the present article, we will mainly restrict ourselves to some particular examples convenient for numerical 
experiments, in which the 2D character of the pump wave incoherence is modeled by the introduction of two parameters 
Awo and Afco, playing the role of the temporal and spatial incoherence, respectively, so that Aj max is similarly written 

as Ai max = CCiAwo + /3iAfco. 

Regarding now under which conditions the pump wave incoherence modifies the parametric coupling, on(^ may first 
easily admit that a sufficient condition for neglecting the pump wave incoherence is that the two following inequalities 
be both satisfied: 

AlTOax <| -Dl I 

^2max <\D2\ {1 1 1. 6) 

In the latter inequalities, Di denotes the dispersion relation 

D^,k. (7, fc) = 7 + ^ - ^kf ) + {1 1 1. la) 
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In the envelope approximation, Di reduces to 

A, k, (7, fc) = 7 - K • Vgi (ki) + lyi {IllJb) 

where, 7 and k denote the time and space growth rates. (It can be noted that our definitions of Di and Di are consistent 
with each other, in the sense that if one considers ab initio the coupled mode equations in their envelope form (11.10), 
Di reduces identically to Di without any further approximation: in this case, one has indeed 0)^- = Ki ■ 'Vgi, from 
which the relation Di = Di follows.) 

Concerning the inequalities (III. 6), we will see more generally that the coherent or incoherent character of the 
parametric couphng is controlled by inequalities between Ai„jax and \ Dj \. On the other hand, since 7 and k are 
related through the dispersion relation, the size of Dj depends upon the nature of the instability which is considered, 
namely convective, absolute, or spatially amplifying. It therefore follows that the coherent or incoherent character 
of the parametric coupling depends itself on the nature of the instability. Defining now as "incoherent" the domain 
where the pump wave incoherence induces a reduction of the parametric growth, one realizes that such a domain 
should be specified as " convectively incoherent" , or " absolutely incoherent" , or " incoherent for spatial amplification" . 
For instance, in the case of a purely temporal incoherence, with Vgo — 00, it has been shown in (Laval, et.al.) that 
the pump wave incoherence reduces the convective growth for AoJq > Max(j/i, V2) and Awq > 70, whereas in the case 
of an absolute instability, the absolute growth can be reduced only if Auq satisfies: 

Awo > 70 I ^ I^^V ' ^ ' ■ 

[For simplicity we consider the case where | Vi \»\ V2 \ and 2^2/ | 16 |>> I'l/ | Vi |.] One thus sees with this simple 
case that the reduction of the absolute growth is much more diSicult to achieve than that of the convective growth. 

We may recall here that the convective growth rate characterizes the instability in the early stage of its development, 
by contrast with its long time evolution which is characterized either by the absolute growth rate (whenever the 
conditions for their existence are satisfied) or by the spatial amplification growth length (in the opposite case). 

In this article we will derive the conditions under which the pump incoherence gives rise to a reduction of the 
convective and of the absolute growth rates and an increase in the spatial amplification length. It will be found that 
generally the conditions for reducing the absolute instabilities are more severe than those for the convective ones. For 
this reason we will concentrate mainly our discussions on the "convectively incoherent" domain, which may thus be 
regarded as the largest domain in which a reduction of the parametric growth may be expected in the early stage 
of the instability development. In order to simplify our terminology, the "incoherent" domain, without any other 
specification, will henceforth refer to the previously defined "convectively incoherent" domain. 

Restricting ourselves to the convective instabilities, one may rewrite inequalities (III. 6) as: 

Aimax < 1^1,10 

A2max< 1^2,70 (///.8) 

in which we used 7 ~ 70 in the coherent regime. As said before, these inequalities represent only a sufficient condition 
for the pump wave incoherence being ignorable; the question is thus whether both of these two inequalities have to be 
satisfied, or only one. In order to answer this question, we first introduce in a brief discussion the validity conditions 
of the statistical equations which describe the time evolution of the quantities < a, > and <| a, |^>. 

B. The Incoherent and RPA domains 

A statistical description for the evolution of the decay waves 1 and 2 is provided by the Random Phase Approxi- 
mation (RPA) equations for the evolution of the average intensities < Oja* >. The latter provide equations for the 

evolution on a slow time - or slow space - scale of the spectral densities n^i^\x,t) of the decay waves [the latter fmictions 

n^*''(x,t)] will be defined in the general case by Eq. (V.2): at this point it is sufficient to say that for a purely temporal 

problem one has the usual relation < a'^\a'^})* >= S(Il — k')nj,'-']. The RPA equations are nonlinear and they may 
account for the pump depletion; the interested reader is referred to standard textbooks[72] for their derivation in 
the general case. In our problem of stability analysis, we face a linear problem; we could naturally use directly the 
RPA equations on which we would make a posteriori the linear approximation which consists in neglecting the pump 
depletion. In this article we will follow a different route. We will first take advantage of the linear character of the 
coupled mode equations II. 6 which will enable us to use the so-called Bourret approximation [70] for the evolution of 
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< fli >; on the other hand, the Bourret approximation for the evolution of the intensities < a^a* > is easily tractable 
within the envelope approximation only, and then only for the special case where the pump wave group velocity VgQ 
is infinite. In the general case it is necessary to make supplementary approximations to the Bourret equations, and 
these approximations are justified in a domain that will be defined subsequently as the "RPA domain." 

Returning now to the condition defining the domain in which the parametric coupling takes place coherently, we 
first define the temporal and spatial growth rates 7<aj>, K<aj>, for (| ai p), and 7(a2>, K(a2) for (| ai p), according to: 

(ai(x, t)} = a^^. (0) exp [ji^^.^t - • x] {I II. 9) 



(I ai 1^ (x, t)) =1 a'^^ p (0) exp 2 [7(^2)^ - ^(^2) • x] 

A factor 2 has been included in the argument of the exponential defining 7<a2> and K<a2> so that the equality 
7<a> = 7<o2> holds in the coherent case. More generally the Schwartz inequality yields the following result 

7(a2> > 7(a0 (7/7.10) 

for i = 1 and 2, in the case of a purely temporal growth (and || K(a2) ||>|| k^^;) II in the case of a spatial growth). 

Admitting that the pump incoherence may only reduce the parametric growth (in the case of a convective instability 
in a homogeneous plasma), one realizes that, whenever either (ai) or (02) behaves coherently, the two quantities 
<| ai p> and <| 02 P> must also behave both coherently. 

In the next section, it will be found that there exists a growth reduction for < 0,1 > only if the mismatch A2 max 
satisfies the condition A2 max > (1^2, 7o)- From the previous considerations, it follows that the wave energies <| p> 
will behave coherently whenever at least one of the two following inequalities is satisfied: 

Ai max < Max(i^i,7o) (777.11a) 



or 



A2 max < Max(z/2,7o) (777.116) 

These incqiialitics define what is referred to henceforth as the " convectively coherent" domain. Accordingly, the 
" convectively incoherent" domain - or simply incoherent domain - is defined by the two conditions 

Ai max > Max(z/i,7o) (777.12a) 

and 

A2 max > Max(z/2, 70) (777.126) 

On the other hand, the validity conditions for the RPA are very stringent. In Section V they will be found to 
consist of two sets of inequalities; the first one corresponds to the inequalities opposite to (III. 6), namely 

Ai max > Max I 7?! I 777.13a) 



A2 ma^ > Max I D2 I (777.136) 

and a second set involving the cross inequalities 

Ai ma. > Max I D2 I (777.13c) 



A2 ma. > Max I Di I (777.13rf) 

Since the RPA dispersion relation yields 7^02^ 7o/min(Aj max) for the convective growth rate, the latter inequal- 
ities reduce essentially to 

Ai max > Max(z/i,7o) (777.14a) 



max 

> Max(i^2,7o) 



(777.146) 
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Ai max > Max(z/2, 7o) {III.Uc) 

A2 max > Max(z/i, 70) {Ill.Ud) 
which can also be recast into the more compact form 

Min(Aj max) > Max(t'i,7o) (7/7.15) 

The latter inequalities (III. 14) or (III. 15), define what is called the (convectivc) RPA domain. 

One can now see that the incoherent domain (III. 12) is subdivided into 1) the RPA domain in which the cross 
inequalities (III. 14b) and (III.14d) are both satisfied, 2) the intermediate domain corresponding to the regime where 
the two inequalities Ai max > (^^1,70) and A2,nax > (2^2,70) are both fulfilled whereas at least one of the cross 
inequalities (III. 14c) and (III. 14d) is not satisfied. 

In this intermediate domain there is no theory which is easily tractable for an explicit computation of the growth 
rate of the average energy <| Ua P>. On the other hand, the Bourret approximation makes it possible to calculate 
the growth rate of the average amplitudes < oi > and < 02 >: it will then be found in Section V that the larger 
growth rate of 7<ai> and 7<02> is given (for the convective instabilities), by the RPA prediction 7^^^ for <| ai p> 
and <| 02 P> within a numerical factor no larger than two. 

For this reason we will argue in Section VI that although the RPA equations are not in their range of applicability 
in the intermediate domain, the convective growth rate 7<a2> for the intensities remains of the order of -y^^^. 
Accordingly, in the whole incoherent domain defined by the two inequalities: 

Ai max > Max(z/i,7o) {Ill.Ua) 

and 

A2 max > Max(i/2, 70) (777.126) 

the convective growth rate 7<a2> can be approximated by the RPA prediction 7^^^. The latter approximation is 
one of the main results of our paper since it makes it possible to eompiite the reduction of the parametric growth 
induced by the pump wave incoherence from the simple calculation of the RPA coupling constants 7^^"* and 7^^^. 



IV. THE BOURRET APPROXIMATION FOR THE AVERAGE AMPLITUDES 

A. Introduction to the Bourret Approximation 

The Bourret approximation[69,70] is a well known equation in the context of propagation in random media; it deals 
with stochastic linear multiplicative equations of the form 

{dt + Lo)A = joSA (IVA) 

where Lq is a linear deterministic operator, S a linear stochastic operator with zero mean, < S >= 0, and A is the 
physical quantity of interest. The Bourret approximation can be simply derived as follows: 

i) by averaging Eq. (IV. 1) one first obtains the exact relation 

{dt +Lo)<A>=jo<SSA> {IV.2) 

with 6A = A- < A>. 

ii) by subtracting the latter from Eq. (IV. 1), one obtains the equation for evolution of the fluctuation 6A 

{dt + Lo) SAo = 70 [.S < ^ > +SSA- <S5A>] 

iii) the Bourret approximation consists then in neglecting the so-called mode coupling term SSA— < SSA > in the 
latter equation; by doing so and neglecting the initial conditions, one obtains for the fluctuation 6 A the relation 
5A = {dt + Lo)~^'yoS < A >, where {dt + Lo)~^ represents the Green's function of the operator {dt + Lq); 
inserting the latter result into Eq. (IV.2), one obtains the Bourret equation for < A>, namely 

{dt + Lo)<A >= 7^ <S{dt + Lo)"' S><A> {IV.3) 

Defining jeS as 7off — 7o ||< S{dt + Lo)^^S >||, it is natural to introduce the effective bandwidth Aweff as 
Aweff = 7o/7eff) and the validity condition for the Bourret approximation can be expressed as 

Alu^s >\ D \ = \\ dt + Lo \\ {IVA) 

The interested reader is referred to Ref. 69 for a detailed discussion of these results. 
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B. The Bourret approximation for the average amplitude < ai > 

1. General three dimensional result. 

By applying the method outhned just above, one easily obtains the equation for evolution of the average amplitude 
< ai >. Consistent with (II 12b), we set 

ai(x, t) = ai (x) exp i (Ki • x — UK^t) + c.c 

and the dispersion relation corresponding to the slow time and space evolution for < ai (x, t) > reads 

(7 + 1^1 -AC- Vgi) =7i {IV.5) 

with 



7i = 71 (Ki) = 7o / .-X n 
J -zA2(ko 



c^konfc„ 



) + (7 + 1^2 - K • Vg2) + e 



{IV.6) 



where 7 and k stand for 7<ai> and K<ai>) according to (III. 9). In deriving the latter equation, we neglected for 
simplicity the slow dependence of Ukoik^ and u^J upon ko and k^. For the sake of clarity, we also restrict ourselves 
here to the case of exact resonance (ki = Ki); the general case ki 7^ Ki will be considered in the next Section in the 
RPA context; (consistent with our other notation, 7^^ denotes the coupling constant 71 evaluated at resonance). 
In deriving Eq. (IV.6) we made the envelope approximation for the slow space dependence only, by using ~ 

oj^. + IK ■ Vgi; on the other hand the computation of the resonance mismatch A2 corresponding to a fast space 
dependence is not restricted to the envelope approximation; lastly, the symbol e in Eq. (IV.6) represents the usual 
prescription for the Laplace transform contour. The dispersion relation corresponding to the evolution of the average 
amplitude < a2{x,t) > can be written in a similar way as 

7 + 1/2 - K • Vg2 = 72 {IV.7) 

with 

J -lAi (ko) + (7 + z^i - K- Vgi) + e 
Defining the effective bandwidths Awj eflf as 

Awi efi(7, «) = 7o/72 {IV.9a) 

and 

Aw2 eff(7,«) = 7o/7i (IVM) 
The validity condition of the Bourret approximation for the average amplitude < ai > reads[16], 

Aw2eff >| £'2 | = |7 + !^2-K-Vg2 I {IV - 10) 

and 

7^ < Awi eff Aa;2 eff {IV - 11) 

the same naturally holds for the evolution of < a2 > with (1 ^ 2). 

2. Markov limit 

At this point let us consider the so-called Markov limit of the coupling constant ; the latter consists in taking 
the limit | £>2 | = | 7 + ^^2 - « • Vg2 |<< A2max, and in using (-iA2(ko) + e)"^ - iPPA^^(ko) + 7n5(A2(ko)). 
In this limit the coupling constant 7^ is given by 

7f = 7o' / rfkon£V5 {62 (ko)) (7^.12) 
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where the superscript M stands for "Markov;" the same expression holds for 7^ with (1 ^ 2). The latter two 
expressions for and 7^ are identical to the RPA coupling constants ^ff"^ and ^22"^ to be derived later. At this 
point it is sufBcient to remark that whenever the Markov limit can be taken, the orders of magnitude of and 7^ 
are 

7f ~ 70VA2 „,a. (/V^.13) 



7f ~ 7oVAi ma. 

These estimates follow simply from the definition (IV-12) for 72^ and from the normalization condition f dkon'j^^ = 

1. On the other hand it appears to be convenient, when discussing the continuity between the coherent and incoherent 
results, to rewrite the inequalities (III-12) limiting the incoherent domain in terms of the spectral widths Awi and 
Au)2 defined by the following relations 

Aa;2 = 7oV7f =7oV7fr^ (IV-U) 



According to the estimates (IV. 13) the order of magnitude of /^uJi is given by 

(IV.15) 

It can also be seen that the quantities Acoi can be expressed as Aui = Atc'jeff(7 = 0, k = 0), so that the validity 
condition (IV-10) for the Bourret approximation for < ai > becomes, in the Markov limit, Auj2 >\ D2 \. The 
spectral widths Aui are given in Section VII for the examples of the pump wave correlation function corresponding 
to numerical solutions. Their general expressions are computed in Ref. [7] in terms of the spectral width parameters 
Afco and A^o for the case of interacting wave- packets; it will be seen that they can be expanded in a similar way as 
Aimax, namely 

Auji = tiAko + CTiAGo + P^A0l {IV.16) 

where the parameters ij, ai and /3i can be expressed in terms of the waves parameters duj"^ /dK.j and d^ui^ / d'K.jd'K.ji . 
The latter parameters are all of the same order of magnitude as the corresponding t^, Ui and /3j of the expression 
(III-2) for Ai maxj although they may differ nonetheless by a numerical constant of order unity resulting from the 
averaging procedure over fco involved in Eq. (III-23), in all cases however the ordering 

AuJi = O (Aj max) 

holds and the quantity Aj max will henceforth be replaced by Aw, in the inequalities defining the domain of validity for 
the incoherent results. Moreover the convective growth rate, well above threshold, is easily found to be 7 ~ 7g / Aw2 so 
that the validity conditions IV-10 and IV-11 of the Bourret approximations for < ai > and < 02 > in their Markovian 
limits become simply 

Aa;2> 1^2,70 {IVAla) 



Aa;i> 1/1,70 (/V.176) 

The latter inequalities justify the conditions (III-12), written in the previous Section in terms of Aui eff for simplicity, 
defining the incoherent domain, that is to say the domain where the Bourret approximation is correct for the two 
average amplitudes < ai > and < 02 >. 

Concerning now the validity condition for the Markovian limit itself, made when deriving the coupling constant 
'yf^ = 7q/Acl;2, there is no general criterion, except in the special case where the envelope approximation is correct - 
in this case the validity condition III-27a for the Bourret approximation for < ai > in the convective regime justifies 
a posteriori the Markov limit Aa;2 » D2- This result is discussed in the next subsection; it is shown in particular 
that in the case where the envelope approximation is valid, there is no intermediate regime between the incoherent 
domain IV-17 and the coherent domain (III-ll); namely, either the Bourret approximation is correct and the coupling 
constants take their Markov limits 7j^ = 7o/At<;2, 7o/Aa;i in the incoherent domain, or the coherent results apply. 
In addition, the two sets of results are continuous from one domain to the other. 

On the other hand in cases where the envelope approximation is not correct, there is the possibility for an inter- 
mediate regime where neither the Markov limit of the Bourret approximation nor the coherent result apply; such a 
situation is discussed in Ref. [7]. 

In order to illustrate the continuity between the coherent and the incoherent domains in the case where the envelope 
approximation is correct, we consider in the next sub-section the special case of a Kubo- Anderson process. 
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C. The Special Case of a Kubo-Anderson Process 

The Kubo-Anderson Process (KAP) is an example of a stochastic process for which the Bourret approximation is 
exact, whatever the spectral width Awq = Vgo^ko is. In the case of the coupled mode equations, it makes it possible 
to compute the coupling constant 7i (and therefore the growth rate) as a continuous function of the spectral widths 
AuiQ. The interested reader is referred to the references[14, 15,69] for an introduction to the Kubo Anderson process. 
In a one-dimensional geometry for the pump wave, and in the limit of the envelope approximation, the spectral density 
of the pump wave, when it is modeled by a KAP, is given by 

(0) ^ n-^Ako 

'° (ko-Kof + Afc2 

Performing the integration over ko in Eq. (IV-6), one obtains = j^j Auji eS with 

Aw2eff = 7+i^2-K-Vg2 + Aw2 (JV".19) 

where Aa;2 is given by 

Aa;2 =1 VgKs - Vgicosax \ Ako = m^^^o {IV.20) 

The quantity ?72 =| 1 — (V^2/V^o)cosai | generalizes the quantity r]2 defined previously in a similar framework. [16] 
In the latter expressions, VgoAko has simply to be replaced by Awq in the case of a pump wave with an infinite group 
velocity Vgo = oo; this limit corresponds to a model in which the pump wave fluctuates in time only, i.e. where 
the function S{x — Vgot) in the envelope equation II-9 is a function S{t) of time characterized by a correlation time 
Tc = Aw^^. Lastly, in order to make a connection with the previously introduced definitions, the expression (IV.20) 
for Ati;2 can be seen to correspond to t2 = *2 =| Vgo — V^2Cosq!i |, following our general expression (IV-16) for A0J2 (in 
the case of a one-dimensional geometry for the pump, one has naturally A^o = 0). The quantity AuJuff is defined in 
a similar way by substituting (1 ^ 2). 

In order to discuss the continuity between the coherent and the incoherent domains, it is convenient to introduce 
as before the quantities Di defined by 



Di=^ + Ui-K-Vgi {IV.21) 



The dispersion relation for < ai > is thus 



where for < a2 > it is 



Di = -^-^ {IV.22a) 

D2 + AiJ2 
2 

D2 = -^r^ (IV.22b) 

Di + Au}i 

As stated before, these two dispersion relations are exact for a KAP and they make it possible to investigate the 
behavior of the average amplitude < > as a function of the spectral width Aluq. The first remarkable point is, as 
first noticed by Thomson and Karush[14], that the two dispersion relations (IV-22) are not identical. Due to phase 
mixing, < > does not necessarily behave like <| |>, nor < ai > like < 02 >. 

For simplicity we now restrict ourselves to the convective instability in our discussion concerning the behavior of 
the average amplitudes < ai > and < 02 >• After setting k = 0, one can check directly, using the exact dispersion 
relation (IV-22a) for < ai >, what has been announced just above. First, the validity condition Aa>2 cff >| D2 \ of 
the Bourret approximation necessarily reduces to the condition Aa;2 >| D2 |, and justifies a posteriori the Markov 
limit 7i = 7^^. Since the latter condition Ati;2 > -D2 | is itself equivalent to Aa;2 > (70,1^2), it is natural to define 
the incoherent domain for < ai > as the domain Auj2 > (70, ^^2)- (The incoherent domain for < a2 > corresponds 
naturally to the domain Aoji > (70,1^1).) Second, there is no intermediate regime between the coherent domain 
{Au)Q — > 0) and the incoherent domain for < oi > in the sense that there is no discontinuity between the usual 
dispersion relation of the coherent case {Di = 7o/£'2) and the one corresponding to the Markov limit of the Bourret 
approximation {D2 = 7o/ I A0J2 |); consequently the growth rate and the threshold are continuous from one domain 
to the other. The same naturally holds for < 02 >. 

What can we infer regarding the behavior of the intensities <| p> from the results concerning the average 
amplitudes? This problem has already been investigated by several authors[14-16]. As explained above, from the 
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inequality 7{a2) > 7(a) j one may assert that the system behaves necessarily in a coherent way in the coherent domain 
Auji < (7o,j/i) or Alo2 < (701^^2)- On the other hand, in the complementary domain defined previously as simply 
"the incoherent domain," one may only use the latter inequality: first, it follows that a lower bound for the threshold 
for the intensities is given by the lowest threshold for the average amplitude, namely 

7^ > Min [Aa;i!^2, Aa;2Z^i] {IV.23) 

secondly the growth rate 7<o2> satisfies the inequality 

7(„2) >Max[7(„,),7(„,)] {IV.2A) 

which, well above threshold, reduces to 

7(a2) > 7o/Min (Acoi, Ac^a) {IV.25) 

in a domain where AcOi > 70 for i = 1 and 2. The question is thus whether these lower bounds are actually attained 
in which case it would be sufficient to simply consider the most unstable average amplitude. In order to answer the 
question we first have to derive statistical equations for the wave intensities, which we undertake in the next Section. 

V. THE RPA EQUATIONS FOR THE WAVE INTENSITIES 

The Bourrct technique can be applied exactly in the limit VgQ = 00 in order to derive statistical equations for the 
wave intensities. In the case of a Kubo Anderson Process these equations are exact for any spectral widths Aloq, so 
that they make it possible to investigate the limit in which they take the form of the RPA equations. For the general 
case Vgo ^ 00 additional approximations have to be made to the Bourret-like equations in order to obtain the RPA 
equations. 

A. The Bourret approximation In the limit Vgo = 00. 

Let us consider the coupled mode equations in their envelope form II-9, in which the stochastic function 5 (x — ^got) 
is a Kubo-Anderson process S{t); the latter is characterized by a correlation function of the form {S{t)S{t + r)) = 
exp (— Awo It]), corresponding to a spectral density in frequency space given by 



.(0) _ ^^0/^ 



AlUq + u;^ 



(V.l) 



In order to derive the Bourret equations for the wave intensities, following Ref. [16], we consider the two point 
correlation function < ai{x,t)aj{x' ,t) > with i and j = 1, 2. The spectral densitites are defined as 



{ai{x,t) 



a*{x,'t)) = J dkjn« exp ik^ • (x- x') {V.2) 
where we allow for a slow time and space variation of the spectral density h^} . In order to be consistent with the 

i 

previous definition (III. 9) of growth rates, we set 

fig (x, t) =1 a^J^^ (0) f exp 2 (7t - K • x) {V.3) 

where 7 and k stand now for 7(0=) and K(^a^y 

The Bourret approximation is easily derived for the set (ai(x, t)a*{x', t)) to give the following system 



= 7fi (k'l) n\:/ + n'X' (^•4«) 



b2h'i) = 72^2 iK) [n\:/ + n'^L,\ {VAb) 
where Di is defined as before as Z), = 7 + — k • Vg^; the coupling constants -y^ are given by 

7g (bi + D2 + Aojo) 

ll (kO = ^-2 (^-5) 

+ ^2 + AWoj + • (Vgl - Vg2)) 

and the superscript stands for "Bourret" . 
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i) Let us first consider the Markov limits of the latter expressions and show that they exactly reproduce the RPA 
results. The Markov limit of corresponds to the domain | Di + D2 \« Aloq, in which case the coupling 
constants 7,f take the following form 



B,M It i\ 



(F.6) 



Aa;2 + (k^(V,i-V,2)f 

On the other hand it will be seen in the next subsection that the RPA equations take the general following form 



D2n^^l = ll J <n^%5 (A^ k^)) 



k' — k' 



k' — k' 



The quantity A2 (kQ,k'j), consistent with our notations for the primed quantities, is given by 

A^ (kj,, k'l) = A2 (Ko + k^,, Ki + k'l) 



{V.7a) 
{V.7b) 

(KB) 



where A2 (ko, ki) has previously been defined as A2 (kg, ki) 
limit A2 (ko,k'i) reduces to 



oJ^fl —(jJi^^' In the envelope approximation 



,(1) 



A^ (k^,, k'l) = k^, • Vgo - K ■ V,i - 04> - k'l) • V,2 



(K9) 



[It may be noted again that these definitions of Aj and Aj are consistent with each other in the sense that, 
if one considers ab initio the coupled mode equations in their envelope form, A2 (ko,ki) reduces identically to 
A2 (kg, k'l) without any further approximation where k^ = kg, — Kq,. The same definition holds for A'^ (kg, kj) 
with the substitution (1 ^ 2). 

In the somewhat degenerate limit Vgo = 00 considered in this subsection, kg • V^o has to be replaced by w'q, 
corresponding to the limit kg = in A-; accordingly, in the integral appearing in the RHS of Eq. (V.7), the 

quantities nj^^dkQ and A ■ (kg, k^) have to be replaced by h^^dujQ and Aj (uj'q, k^) respectively where A^ (wqi kj) 

denotes A^ (kg, k^) in the limit Vgo ^ 00, kg with Wg = kg • V^g. One easily finds 



A2 (4,k;)=4-k;.(v,i-v 



S2j 



(KIO) 



and (1 ^ 2) for A^ 



It is also only in the same limit Vgo 00 that the integration over kg in Eq. 



(V-7) does not involve the spectral densities n\^i and ny^^_^^, ; the RPA equations take then a simpler form, 
namely 



I>rn? 



7iT^ (kl) 



,(2) 



where the coupling constants are given by 



(1) , (2) 



{V.U) 



722 



^^^(k'2) 



Performing the integration over wg in the latter expression, one finds 



7 



RPA 



(K12) 



(K13) 



The RPA equations correspond thus exactly to the Markov limit of the Bourret approximation for the intensities 
(lap). 
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(ii) Wc now consider the dispersion relation corresponding to the exact Bourret equations (V-4). Taking k2 = — 
in Eq. (V-4b), one obtains 

D1D2 = 7n (k'l) (^1 + ^2) (V.14) 

The latter dispersion relation generalizes to three dimensions the results obtained previously by Laval et al.^^ 
Following their analysis, one may remark that the coupling constants 7^ (k') = (^k^) reach their maximum 
7fj = 7^2 for k'^ = 0, corresponding to the exact resonance conditions ki = Ki and k2 = K2 for the two wave 
packets considered here. In this case the dispersion relation (V-14) reads simply 

b,b,=^l . (K15) 

from which it follows that the dispersion relation -D1-D2 = 7o of the coherent case is recovered for Awq <| 
D1+D2 I- The opposite limit Awo >| D1+D2 | corresponds to the validity domain for the Markov approximation, 
and therefore for the RPA equations. There is thus no intermediate domain between the coherent and the 
RPA ones in the case Vgo — 00. This result can be tracked back to the fact that for Vgo = 00, one has 
Auji = Auj2 — AujQ, as can be seen from the general definition (lV-20) for Acj^, so that the cross inequities 
(III-14c) and III-14d) are automatically satisfied in the incoherent domain (III-12), and consequently the RPA 
equations are exactly applicable in the whole incoherent domain. To summarize, the parametric coupling takes 
place coherently in the coherent domain 

Awo < Max(i/i,7o) {V.16) 
whereas in the complementary - or incoherent - domain 

Awo > J^i, 1^2,70 {V.17) 
the coupling constants jfi , 7^ may be reduced to their RPA values 

#^=7r^ = 7f = ^ {V.18) 
for i = 1 and 2. The corresponding dispersion relation is then simply 



iii) Regarding now the comparison between the average amplitudes and the intensities, one may easily obtain the 
threshold and the growth rate for the intensities from the exact dispersion relation III- 14. 

The threshold reads 7g > z^i ^2(1^1 + 1^2 + Aa;o)/(;^i + 1^2)', it reduces to the coherent result 7q > 1^11^2 for 
Au)o < vi + 1^2 — Max{u\,y2) and to the RPA result, namely 7q/A(jJo > J^i 1^2/(2^1 + ^2) ~ Min (i^i, ^2) in the 
opposite case. These expressions can be easily checked to be identical to the ones obtained from consideration of 
the average amplitudes, i.e. by looking at the most unstable average amplitude, as explained in Section III-2-2. 

Concerning now the growth rates, the growth rate for the intensities is given, well above threshold, by the usual 

coherent expression 7<a2> = 70 for 70 > Ati;o and by the RPA one 7<a2> = 27q/Acjo in the opposite case. 
It is interesting to notice that the correct value for the growth rate (27<a2>) of the intensities is exactly twice 
the lower bound (IV-25) obtained from consideration of the average amplitudes. The square amplitudes behave 
indeed as (ai)^ ~ («2)^ ~ exp (2 gammaQ/ Aujo)t, whereas the RPA equations predict in the same incoherent 
regime (| ai p) ~ (| 02 P) ~ exp(47Q/Awo)t. This result is most easily understood by remembering that, in a 
stochastic process, such as this one, no two realizations are exactly the same, and there is probability associated 
with a given mode's history in amplitude space. Given that probability distribution function, one can compute 
the average mode amplitude or the average mode intensity. If the distribution is narrow, then the intensity 
growth rate will be close to twice the amplitude growth rate. However, in this case, the fluctuations about the 
mean are comparable to the mean in the incoherent limit so that higher order moments of the mode amplitude 
will have a more than proportionately larger growth rate. This facet of the problem is discussed more completely 
in Sec. VII. This resiilt is not unexpected and is actually a general feature of the fully incoherent regimes: in 
an incoherent regime, the growth rate of an average amplitude is indeed given quite generally by 7<a> = 7o''c) 
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where 70 is the coherent growth rate and Tc the correlation time of the stochastic process. Since the coherent 
growth rate of the quantity | a p is 270, we deduce that the growth rate of the average intensity is given by 
7(^2) = (270)^ Tc = 2 (27^„)) , corresponding thus to twice the lower bound '^"y{a) obtained from Eq. III. 10. More 
generally it will be argued in Section VI that the thresholds and convcctive growth rate regarding the intensities 
can be obtained quite accurately from the results corresponding to the average amplitudes, whenever one has 
711 >> 722 or 722 >> 711, and at worst within a factor of two in the case 711 ~ 722- 

B. The RPA equations 

In this subsection we derive the RPA equations in the general case from the Bourret equations for the intensities. 
Following the definition (V-2) for the spectral densities, we set 

(ai (x, t) al (x', t)) = j dkin]^l ("~^~~' ^'^^ ' ~ ^'^ 

and we Laplace transform in space the slowly varying spectral density n^^ (x, t) by setting n^^ (x, t) = J dn exp — 
2k ■ X n^^{K,t). The quantity n'j^^{K,t) is easily checked to be given by 

The equation for the evolution of a)^^_^_^^ ('^ki-ire) obtained from the coupled mode equation (II. 7a) and 



reads 



70 j dkoa£^ (ak„_(ki+i«))* 



+ [c ■ C, K ^ —k] 



{V.21) 



where we neglected as before the slow variations of i/^j and I'^g^ki ko-ki upon ko and k. Similar equations can be 
written for the quantities 



"ko-ki+iK I "ki-iK J 



^k2+iK l"k2-iK 



so that the complete set for the quantities aiaj* has the form of a stochastic linear multiplicative system (where the 



random variable is the set a"^^), upon which the Bourret approximation can be performed. It is convenient at this 
stage to Laplace transform in time the resulting equation by setting 

»L'.'(^.')-/'i7»L'.'(«,7)exp2., 

and similarly for the other quantities. The Bourret equations take the form 



(ikon^"'i?2 (ko,ki) 



where the resonance function R2 (ko, ki) is given by 

ii2(ko,ki) = 



"k/ ('^' + "ko-ki 7) 



DI + D2 



{Di + D2f + iA2 (ko,ki))' 



(V.22) 



(y.23) 



and where the operator Df{'y, k) denotes the Hermitian part of Di (7, k), i.e. 



D+ (7,k) = -[Z}i (7,«) + (Z?i (7*, «*))*]. 
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In the general case is given by 



(7, 



7 + 



i (^u 



ki -\-iK, 



,(1) 



+ 1^1 



(V.24) 



which reduces to D^]^^ = -Di,fei = -y + fi — k ■ Vgi (ki) in the envelope approximation. As a matter of fact, the 
envelope approximation can be made for the slow variation associated with the variables 7 and k; thus the Di and 
D2 appearing in the definition of R2 (ko,ki) can be reduced to -Di^ki and and £'2,ko-ki respectively, and similarly 
the operator can be replaced by -Di,ki (by contrast the quantity A2 (ko,ki) contains the fast dependence on the 
integration variable ko and cannot be reduced to its envelope approximation in the general case.] 

Making the envelope approximation for the slow dependence, one obtains the following set of equations. 



(^i-7fi(ki) 



70^ 



dkoi?2(ko,ki)n|,°>['^' 



(£)2 - 72^2(k2)) n^t^ = ^2 1 ^ko^, (ko, ki) n«n£,)_ 
where the "diagonal" coupling constants 7,^ are given by 

7fi (fci) = 7o / dkong^^2 (ko, ki) 

with 

^2(ko,ki) = 



ki 



D1+D2 



(D1+D2) +(A2(ko,ki))' 



(V.25) 



(V.26) 



(V.27) 



and (1 ^ 2) for 7f2(k2). 

The set (V.25) constitutes the Bourret equations for the intensities (| ai p) and (| 02 P). This set is still very 
difficult to solve since, in addition to the convolution integral on the RHS, it contains various quantities which are 
non-Markovian. The coupling constants 7^^ can be further simplified in the Markov approximation, denoted by 'y^'^ ■ 
In this approximation, i.e. | A2 |>>| Di + D2 \,R2 {ko,ki) can be replaced by nS (A2(ko,ki)) so that the system 
(V.25) becomes 



= Jo I dkon^°Ks (A2(ko, ki)) [nW + n^tk. 



{V.28a) 



*k2 



70^ 



dko4°>^(Ai(ko,k2)) 



,(2) 



'ko-k2 



{V.28b) 



The latter set constitutes the RPA equations for the coupled mode equations neglecting pump depletion and written 
in the Laplace variables 7 and k. In the real space the operators Da have simply to be replaced by 

I {dt + Vga ■ Sx + 2!/„) . 

In order to make a closer connection with the special case considered in the previous subsection (V.7), one can 
easily check that the latter set reduced identically to the system (V.7), when dealing with primed quantities (i.e. after 
factorization of the fast variation). One may also remark that the mismatch A2 (ko,ki) reduces to A2 (ko,ki) = 
ko • Vgo — ki ■ Vgi — {ko — fci) • Vg2 in the envelope approximation, in which limit the set (V.28) can be recognized to 
be the 3-D generalization of the Laval et al. results. "'^^ 

The validity condition for taking the Markov limit of the Bourret equations (V.25) is Min Awj >\ Di + D2 \ — 
Max; I |. On the other hand, by using the fact that well above threshold the growth rate is of the order of 
7 ~ 7o/Min(Aa;i), one easily finds that the latter inequality reduces to 



Min(Aa;j) > Max(z/i,7o) 



(V.29) 
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that is to say the condition (III. 15) previously written in terms of AjMax- This condition entails in turn the inequality 
7q < AwiAw2 which ensures the validity of the Bourret approximation. Therefore, the condition (V.29) defines 
entirely the RPA domain. 

The RPA set (V.28) remains a system difficult to solve due to the integration over ko on the RHS. One can however 
make some additional approximations which are consistent with the validity condition (V.29). By assumption the 

spectral density n^^^ of the pump wave peaks when ko is equal to the mean wave number of the pump wave Kq; 
consequently, for a for a given angle of observation defining the direction of the scattered wave, the RPA coupling 

constant 

7ir^(ki) ^ 7n'^(ki) = 7o' / dkon^^nS (A2(ko, ki)) (y.30) 
is maximum for ki = Ki, corresponding to the exact resonance. In this case 7/l^'^is given by 

^RPA ^ ^fiPA(K^) ^ ^2 I rfkon(%5 (A2(ko)) (K31) 

corresponding to the result announced in Section IV, namely 7^^*^ = 7^. If one now takes ki = Ki in Eq. (V.28a) 

the relation A2 (ko, Ki) = defines a surface in ko space that contains the point ko = Ko(inl — D, A2 (ko, Ki) = 
would require ko = Kq uniquely; as a matter of fact, the function n^^ attains its maximum at this point. Due to the 

large spectral width imposed by the RPA validity conditions (V.29), the function nl^^_j^^ is a slowly varying function 
of ko, so that one does not make a significant error in the integral appearing in the RHS of Eq. (V.28a) by replacing 

'^k^o-Ki ^'•'^^'^^ '^k! approximation which is actually exact in ID), and similarily n^g_j<^^ by in Eqs. (V28b), 
and obtaining 

£>ini =7fi^^ (ni+na) 

^2^2 =7^2^^ 712) (1^.32) 

where the RPA coupling constants 7^^^ are given by 

7^^ = 7oVAa;i 

and 

72Y^=7oVAa;2 

according to our definition (IV. 14). The quantities n„ denote, for the for the sake of simplicity in the notations, the 
spectral density n'^^ at the point of exact resonance k^ = K^. Equations (V.32) provide the desired system with 
which the stability of the intensities <| p> is analyzed in the remaining of the paper. 



VI. THRESHOLDS AND GROWTH RATES FOR THE CONVECTIVE AND THE ABSOLUTE INSTABILITIES AND 
SPATIAL AMPLIFICATION 

In this section we derive the expressions of the convective and absolute growth rates for the parametric instabilities 
as a function of the two quantities Aa.-'i and Acj2 characterizing the pump incoherence; we also compute the rate of 
spatial amplification. Let us first recall the physical meaning of these various growth rates by considering the Green 
function, that is to say the response of the parametric coupling to an infinitesimal impulse given to the system at t 
— and X = 0. For simplicity wo restrict ourselves to a one dimensional problem, and we denote as Xm{t) the point 
where the Green function G(x,t) is maximum. 

For an infinite system, the convective growth rate, denoted as 7^^°""^, is defined to be the time growth rate of 
Goo{xm{t),t), i.e. 

Goo(a;™(t),t)r-.ooexp(7'°""i) 

Quite generally the point Xm{t) moves in time according to Xm{t) = Vmt, where V^(= Fi + V2)/2) is the group velocity 

of the maximum. For a finite system, the convective growth rate characterizes therefore the transient behavior only; 
although the Green function Gl for a plasma slab of length L does not reduce to the Green function Gb of the infinite 
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case, it is usually admitted that the infinite system convective growth rate 7^°""^ does properly describe the early time 

behavior for t< LjVm- 

Regarding now the long time behavior of the parametric coupling, one has first to look for the existence of absolute 
instabilities. For an infinite system, when such instabilities exist, the absolute growth rate is the time growth rate of 
Goo{x,t), i.e. 

\G^{x,t) I Jl^^exp - ^-''^a;] 

for finite values of .t; in the latter expression k"'"* denotes the space growth rate associated with the absolute instability. 
For a finite system, there exists in general a critical length Labs above which there are unstable normal modes, so that 
the system will behave asymptotically in time as the most unstable normal mode. Similarly it is usually admitted 
that whenever the plasma length L satisfies the condition L » Labs, the growth rate of the most unstable normal 
mode is well approximated by the absolute growth rate 7"''* of the infinite case. 

When absolute instabilities do not exist but the system is convectively unstable, the long term behavior is determined 
by spatial amplification. The latter corresponds to setting a source Soexp — iust at x = and looking for the tOg 
which maximizes the spatial amplification rate k{cOs)', Bers and B riggs have indeed shown that in such a case the 
time asymptotic behavior of an infinite system in the presence of the source goes like 

[GooSoexp {-iiJst)]^^^^^ ~ exp {eK{uis)x) 

where e = ±1, depending upon the direction for spatial amplification. Such behavior also indicates that a significant 
spatial amplification may be expected in a finite plasma slab whenever the condition | k^'^L |>> 1 is fulfilled, where 
k"'* denotes the maximum k{u)s). 

Lastly there is also the possibility for a mixed situation where there is the coexistence of an absolute instability and 
spatial amplification. Indeed such a case will be seen to be provided by the R.P.A. dispersion relation. In this case the 
short and long time behaviors are still given by the convective and absolute growth rates, for the reasons explained 
just above. In the case however where the spatial amplification growth rate | k^" | is significantly larger than the 
space growth rate | Kabs \ associated with the absolute instability, one may expect the existence of an intermediate 
regime in the time evolution during which the spatial amplification feature dominates that of the absolute instability. 



A. Convective instabilities. 

1. Domain of applicability of the RPA equations. 

To begin with, let us first justify our conjecture concerning the applicability of the RPA results in the intermediate 
domain where neither the coherent nor the RPA equations could a priori be used. To do so we compare the RPA 
predictions with those for the average amplitudes. 

The dispersion relation corresponding to the RPA equations V.32. is 

(j^i - 71) (-D2 - 72) = 7172 {vi-'^) 
and the dispersion relation for the average amplitude (oa) Eq. (IV. 5) is 

£>« - 7a = {VI. 2) 

where 7^ denotes the RPA coupling constant ^^a"^ exact resonance, namely 71 = 7Q/Aa;2 and 72 = 7o/Ati;i. The 
spectral widths Aloj are defined by Eq. (I-IO) and (IV. 20). 

As explained in Section III, the intermediate domain is defined as the domain where equations (VI. 2) are correct 
for the two waves a = 1 and 2 (i.e. the diagonal inequalities Aua > Ma,x{i'a,jo) are both satisfied) whereas the RPA 
equations are a priori not applicable (because one of the cross inequalities Acua > ^/s, with a 7^ /?, is not satisfied). 
In this domain one can only use the inequality (III. 10), namely Max(7(Q.-)) < J{a^) which states that the actual 
growth rate 7^02^ is no less than the most unstable average amplitude (the latter being reduced as compared with the 
coherent growth rate since Eq. (VI. 2) predicts a reduction of the parametric growth in its domain of applicability). 
The question is thus whether this lower bound is attained or whether 7(a2) remains of the order of the coherent growth 
rate. To answer this question, it is natural to first compare the RPA predictions with those for the average amplitudes 
in their common domain of applicability, namely the RPA domain. It will be seen in the next subsection that they are 
identical with regard to the convective threshold and convective growth rate to within a numerical factor no larger 
than two. 
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The reason for this identity can be tracked back as follows. The RPA dispersion relation VI. 1 can be approximated 
by I?i — 7i ~ if 72 << 71 or by 1)2 ^ 72 — if 72 >> 71, so that the two quantities (| Oi p) ~ (| 02 P) have 
necessarily a growth rate 7(^2) which is of the order of Max 7(ai)- More precisely, denoting by 7_rpa the growth 
rate obtained from the dispersion relation VI. 1, and by 7max the largest convective growth rate given by the Bourret 
dispersion relation VI. 2, it is the matter of a simple calculation to show that the following inequalities 

7max < IRPA < 27 max 

hold formally, i.e. independently of the applicability of the RPA and Bourret equations (the upper inequality occurs 
in the case 71 = 72). Since in the RPA domain one has 7^02^ = ^rpa and Max 7(„.^ = 7max, one deduces that the 
latter inequality can be written 

Max (7(„,)) < 7(„2) < 2Max (7(„,)) {VIA) 

in the RPA domain. This result means that the phase mixing effects on (oq,) arc negligible, at least for the most 
unstable average amplitude. Since this property is satisfied in the RPA domain, which can be regarded as the most 
incoherent domain, i.e. the domain in which the phase mixing effects are expected to be the more important, we 
make the conjecture that the inequality (VI. 4) is satisfied in the intermediate domain as well. On the other hand, 
in this domain one has Max (7(0^)) = 7max; so that this conjecture means 7max < 7(o2) < 27max; on comparing this 
inequality with (IV. 3), one realizes that one can approximate in the intermediate domain, the actual growth rate 
7(a2> by the RPA prediction -)rpa-, with an error no larger than a factor of two. As anticipated in Section III, these 
arguments allow us to apply the RPA results in the entire incoherent domain Awi > I'l , 70 and Aa;2 > J^2) 7o including 
the intermediate domain. 



2. Convective instabilities expressions 

The threshold 7oconv for the convective instabilities is easily obtained from the RPA dispersion relation (VI. 1); it 
reads 

1 /2 

7oconv = [vii'2 (z^i/Awi +zv2/Aa;2)] (^-'"•5) 

and can be approximated by 

7oconv = [Min (Aa;i!/2, ^^2^1)]^'^ {VI. &) 

It is easily checked that the latter approximation corresponds to what can be deduced from the behavior of the 
average amplitudes, [namely from the condition 7max > 0]. We may also notice that the expression VI. 6 for the 
threshold in the incoherent domain goes continuously into the usual coherent threshold expression 70 = i'iV2 at the 
boundary between the two domains, namely Awi = vi and Aa;2 > i'2, or Aa;2 = 1^2 and Awi > vi (again within a 
numerical factor no larger than two). 

The convective growth rate, well above threshold, is given by 

*T=^2(ALr+i) <^"' 

which is of the order of 7inax = 7o/ (MinActij), again within a factor no larger than two [one has 7(o2> = 27max for 
AoJi = Au)2 only, which occurs for instance in the case V^o = 0° encountered in Section III. 2]. One can again show 
that the expression for the growth rate in the incoherent domain goes continuously into the coherent one 7^02^ = 70 
at the boundary Awi = 70, and Aa;2 > 7o and Awi > 70 (once more within a numerical factor no larger than two). 



B. Absolute instabilities 

The absolute growth rate is obtained from the RPA dispersion relation VI. 1 by looking for a double root for k 
and by imposing the well-known Bers and Briggs criterion, i.e. tracing the roots in the K-plane to insure that the 
double root arises from roots coming from opposite sides of the real K axis as one varies the growth rate from 7^^!'^ 

to 7^^2^^ + 00. One finds that this occurs for V1V2 < and that the absolute growth rate is given by 
7^1) = {\V,\ + \V, \)-' 



2 I \2 
7o - (70abs) 



Vl 

Awi 



1/2 



+ 



V2 

AW2 



1/2 



(y/.8a) 
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where 7oabs denotes the absolute threshold 



7o^abs =1 V,V, I ^!^:' '^-'^ {VI.9a) 



Well above threshold the absolute growth rate can be approximated by 



^2 

^<°'> ~ (Max| |)(Min(Aw,-/ \Vj\)) (VIM) 
and the absolute threshold is of the order of 

7o\bs ^1 V,V2 I (Max^) ■ ^^^-^^^ 

Concerning now the domain of validity of the RPA expressions for the absolute instabilities, let us first consider the 

case Alui = Acl'2 = Acjq; in this case the cross inequalities (III.13c-d) reduce to Acjq >| Da \, so that it would appear 
that the RPA validity domain is given by the inequalities Aujq >| -Da | with a = 1,2, where = jf^^'^— n'^^^Vga + i^a is 
computed for 7 = 7^^!^^ (with k evaluated at its double root Kabs)- Actually a detailed calculation [46] performed with 
the exact dispersion relation (V.15) for the case Vgo = 00, shows that the limiting RPA expressions are applicable in a 
smaller region, obtained by imposing the RPA validity condition not at the single point (7abs, '^abs) but on the entire 
path needed to apply the Bers and Briggs procedure: in fact one requires that Awq >| Da | = | 7 — K±{'j)Vga + Va \ 
for 7 varying from the absolute growth rate 7^^2^^ to the value at which one of the two roots K±{'y) crosses the real 

axis. We assume such a procedure is also required in the the general case Ao^i 7^ Aa;2- One then finds | Di |— 7o/ | 
V2 I Min(A(x;i/14) and | D2 |— 7o/ | V\ \ Min(Aa;i/ | Vi |), so that the condition | Aw^ |>| Da \ reduces to the simple 
inequality 

\\Vi\) 

On the other hand, by looking at the behavior of the average amplitudes (a^), one easily finds that in the comple- 
mentary domain Min ( Awi / \Vi\) < 70 / \JV1V2, the parametric system behaves coherently for the absolute instabilities 
and the growth rate reduces to the coherent limit given by 

^abs 2a/| V1V2 I I ^,coh\ (VT^^n\ 

= I y I _^ I I (To - 70absj : {V 1 .Ua) 




(^7.116) 



For this reason we can define the incoherent domain for the absolute instabilities by the inequality (VI-10). 

With regard to the two cross inequalities | Awi |>| D2 \ and AC1J2 >| -Di |, it can be seen that, because they arc not 
necessarily fulfilled when the condition VI. 10 is satisfied, there is again the possibility for an intermediate domain in 
which neither the RPA nor the coherent results apply for the absolute instabilities. One can verify however that the 
single validity condition VI-10 ensures the continuity (within a factor 2) of the RPA expression (VI-8) for the absolute 
growth rate with the usual coherent expression 7^'^'' ~ 270 (Min | Vi \ /Max | Vi l)^^"^ ■ For this reason we conjecture 
that the RPA expression (VI.8) is actually applicable in the whole incoherent domain. As already noticed[16], it is 
interesting to remark that the behavior of (a) and (| a are significantly different as far as absolute instabilities are 
concerned: in the incoherent domain (VI. 10) the average amplitudes are essentially stable, whereas the intensities are 
absolutely unstable as long as 70 exceeds the threshold (VI.9). 

In conclusion, therefore, the absolute growth rate is given by the minimum of 7^^!') given by Eq. (VI. 8a) and the 

coherent growth rate 7^^^ given by Eq. (VI. 11). 
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C. Spatial Amplification 

The spatial amplification growth rate k^'^ is obtained from the RPA dispersion by taking j — iui, with uj real, and 
by looking for a maximum of | Re k | which satisfies the Bers and Briggs criterion [75] for spatial amplification namely, 
the root k in the K-plane has to cross the imaginary axis as one varies 7 from iu to + 00. Although the procedure 
is straightforward, the general expressions of the different unstable roots are tedious due to the large number of the 
independent parameters V^,fa, Aw^. For the sake of simplicity we will restrict ourselves in this subsection to the 
cases where the following ordering is satisfied: 

I Vi |»| V2 I {Vl.Ua) 

»T^ ivi.m) 



\V2\ I Fl 

rer » mi 

Such an ordering is always satisfied in the case of SRS and SBS instabilities and quite often for two plasmon decay. 
It is easy to find that in the above limit there are two roots corresponding to spatial amplification. The spatial 
amplification growth rate for the first root is given by 

= (^7.13) 

and its domain of existence corresponds to the inequalities 

Aa;iAa;2 > 7o (Vl.Ua) 

Awi > ui {VI.lAb) 

Au}2 > 1^2 (VIAAc) 

Auj2i^i < 7o < Awii/2 {Vl.lb) 

The inequalities (VI. 14) are the RPA validity conditions, whereas the inequalities (VI. 15) represent the condition 
for satisfying the Bers and Briggs criterion, i.e., the condition for spatial amplification (this condition implies in 
particular the inequality Aa;2i'i < Aa;ii/2). 

The spatial amplification growth rate for the second root is given by 

= -ji/Au;iV2 {VI.W) 
and its domain of existence corresponds to the following inequalities: 

Aui > vi (Vl.Ua) 

Aw2 > 1^2 {VI. 17b) 

Aa;i>7o|^r/' {VI.Uc) 

V2 

7^ > A(JiU2 {VI.18) 

again where inequalities (VI. 17) are the RPA validity conditions for this root, and as the inequality (VI.18) is the 

Bers and Briggs condition for spatial amplification. It is interesting to remark that the inequalities (VI. 17c) and 
(VI.18), reduce exactly to the conditions (VI. 10) and (VI. 9b) for the existence of absolute instabilities in the RPA 
regime, in the case V1V2 < and where the ordering (VI. 12) holds. In this case one is thus in a mixed situation 
where an absolute instability and spatial amplification coexist together. It can also be easily shown that the ordering 
(VI. 12) implies the inequality | |>| k^^^ |, so that spatial amplification may dominate the absolute instability for 
an intermediate stage in time or for an unbounded plasma. 

Lastly, it is worth mentioning that the two roots and can be interpreted a posteriori as corresponding to 
spatial amplification for the average amplitudes (ai) and (02) respectively. For instance the expression for follows 
simply from D2 = 72, the inequalities (VI. 17) are simply the Markov condition | Da \< Auja, and the condition 
(VI.18) is the condition for instability for the average amplitude (02). We may therefore conclude that the main 
features corresponding to spatial amplification - and to the convective growth rate - for the intensities (| a |^) may be 
obtained from consideration of the average amplitudes (a^), in contrast with the absolute instabilities the features of 
which can be found from the dispersion relation for (| a p) only. 
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D. Comparison of Stability Domains in the Case Vgo = oo 

In this subsection we repeat the various stabiHty domains obtained in the previous subsections in the special case 
Vgo =00, for which one has Awi = Au)2 = Awq, in order to exemplify the different kinds of stabilization that can 
be expected from the pump wave incoherence. To do so, it is convenient to introduce the following dimensionless 
quantities 



Awo = Awo/lo 

V2 = V2/V1 

Wc assume the ordering (VI. 12) to hold, and in addition we suppose ui < U2 so that we restrict ourselves to the 
situations where the two inequalities 

I t;2 |« 1 {VI. 19a) 

Vi < V2 {vi.m) 

are satisfied. 

1. Convective Instability 

In Fig. 2 we reproduce the stability diagram corresponding to the reduction of the convective instability growth 
rate induced by the pump wave incoherence. This diagram has to be understood as representing the effect of pump 
bandwidth Awq as a function of V2, for 70 and v\ < V2 given. As explained in the introduction to this Section, the 
convective instability represents the early or transient time behavior of the parametric coupling. 

The domain of applicability of the statistical theory corresponds to the inequalities 

Awo > 70, 
Awo > Max (yx,V2) = 1^2, 

that is to say, 

AU^ > 1, {VI.20a) 



AuJ^ > V2. {VI. 20b) 

The domain corresponding to these inequalitites is referred to as "VR" in Fig. 2, and stands for "convective growth 
rate reduction;" the solid line presents the boundary between the coherent and incoherent domains. The threshold 
VI. 5 for convective instabilities is approximately 7oconv = (Awqz^i)^^^, so that the pump wave incoherence induces a 
complete stabilization of the convective growth rate in the domain referred to as "VS" in Fig. 2, 

AZJo > z^r^ {VI.21) 

The boundary of the convectively stable region is the dotted line as obtained from Eq. VI. 5. Lastly the parametric 
system is coherently stable in the domain F2 > (corresponding to the usual coherent condition for stability 
7q < Vii'2), the boundary of which is represented by the dashed-dot-dot line. 



28 



2. Absolute Instability 

In Fig. 3, we reproduce the stabihty diagram corresponding to the reduction of the absolute instability growth rate 
due to the pump incoherence, i.e. the growth reduction of the long time behavior of the parametric coupling in the 
case of existence of absolute instabilities. The domain of applicability of the statistical theory of absolute instabilities 
(VI. 10) reads here Awq > 72 | VxlVi 1^^^; on the other hand the threshold (VI.9b) is 7oabs = (AwoJ^2)^^^ The 
intersection between the two lines Aa;o = 7o (^1/^2)^^^ and Awq = 7o/^2 corresponds to = 70 (^2/^1)^^^, which is 

1 /2 

smaller by a factor two than the usual coherent threshold for absolute instabilities, namely v-2. > 270 (V2/V1) ' . Since 
the applicability condition for statistical theory has been obtained from strong inequalities (whereas the threshold 
condition 7oabs = (Aa;ot'2)^^^ results from solving the RPA dispersion relation) the two limits could be reasonably 

1 /2 

made continuous between each other by writing the RPA validity condition as Auiq > 70 {V1/V2) /2 for 1^2 < 

1 /2 

270 (T^2/T4) ' . By doing so, the domain corresponding to a reduction of the absolute growth rate (henceforth referred 
to as the "AR" domain) would be approximated by the inequality 

Acj5J> (2 I V2 1^/^) ^ iVI.22a) 

and the domain for complete stabilization of the absolute growth (referred to as the "AS" domain) would be given 
by: 

AZZJJ^ >| 772 (VI.2'Sa) 

Actually a more detailed calculation [46] based upon the exact dispersion relation (VI. 15) shows that the AR domain 
corresponds to the inequality 

>\ V2 r^/^ (^7.226) 

for F2 <| V2 1^^^, and that the complete stabilization occurs when inequality (VI. 23a) is satisfied in the range 
z^2 <| ^^2 r^^, and when the inequality 

Awo > 2 I t;2 1"^/^ -V2 I V2 

is fulfilled in the range | V2 \~^< V2 < 2 \ V2 \^^'^ . 

In Fig. 3 the boundary of the AR domain is represented by a short-dashed line, and the boundary of the AS 
domain, computed from the incoherent threshold VI. 9a, by a long-dashed line. One may see that the two boundaries 
almost coalesce at V2 =| U2 |^^^, Auja =| V2 |~^^^; this simply means that the AR domain shrinks to zero in the limit 
I W2 1^ for 172 >| W2 1^/^ (for this reason the short-dashed line, the long-dashed fine, and the verticle dashed- dot-dot 
line bound a very small region for T72 >| V2 |^^^)- The verticle line is the absolute coherent threshold Eq. VI. lib. 

It can be easily checked that the domains corresponding to the approximate expression (VI. 22a) and (VI. 23a) do 
not significantly differ from the exact domains (VI. 22b) and (VI.23b), in the range 172 <<| V2 1^^^, i.e. far from the 
absolute threshold. This result justifies in particular all the approximate expressions obtained in the remainder of our 
paper. On the other hand, in the range close to the absolute threshold, namely for | |^^^< F2 < 2 | U2 |^/^, this exact 
solution exhibits the additional AS domain 270 | V1/V2 \^^^ —V2 \ V1/V2 \< Awq < 70/^^2, corresponding to a regime 
where a small Acjq may reduce the absolute growth rate. Such an effect is not really surprising since it means that a 
small frequency bandwidth may have a large effect when the parametric coupling is close to the coherent threshold, 
namely it converts an absolute instability into a convective one (which in turn gives rise to a spatial amplification, 
as it will be seen in VI. D3), by decreasing the usual coherent absolute threshold 772 = 2 | t;2 |^^^ into the absolute 
threshold (VI. 23b) modified by incoherency effects, namely V2 = 2 \ V2 \^^'^ — \ V2 \ Aujq. 

3. Spatial Amplification 

In Fig. 4, we reproduce the stability diagram corresponding to the reduction of the spatial amplification growth 
rate k*". The first root (VI. 13) for the spatial amplification growth rate is given in terms of dimensionless quantities 
as 

M = - (AZJo)"^ (^7.24) 

where the dimensionless growth rates 7« are defined as k = k^°'Vi/^q. The domain of existence of this root, referred 
to in Fig. 4 as "ARl" (for "spatial amplification reduction"), is defined by inequalities (VI. 14) and (VI. 15), namely 



Awo > 1 



(F7.25a) 
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Auo > V2 {VI.25b) 

Auo > F2 ^ {VI.25c) 

The boundary of this domain is represented by the dotted-dashed Une in Fig. 4. 
The threshold condition (VI. 15) can be written as 

Awo < V^-'^ {VI.26) 

The boundary of this domain "AS" is given more exactly by the dotted line as in Fig. 2. 
The spatial amplification growth rate for the second root is 

K2 = - (A^0V2)"^ {VI. 27) 

and its domain of existence, referred to in Fig. 4 as "AR2," corresponds to the inequalities (VI. 17) and (VI. 18), i.e. 

Awo < i^2~^ {VI.28a) 

Awo >| V2 {VI. 28b) 

in addition to inequalities in (VI. 25a) and (VI. 25b). The boundary of this domain is represented by the decreasing 
dotted-dashed line and the dashed line in Fig. 4. It is interesting to remark that in the domain 1/2 > 2 | 1^^^, the 
coherent spatial amplification growth rate is Kcoh = so that k varies continuously from its coherent value to 

Ki = — l/Awo, along the boundary AtUq = V2 (the dotted-dashed line for F > 1), whereas k is strongly reduced from 
Kcoh to Ki along the boundary Acoq = I/F2 (the dotted-dashed line for F2 < 1). 

Similarly, in the domain V2 « 2 \ V2 |^^^, and V1V2 > 0,k varies continuously from its coherent expression 
K = ~l/v2^'^ to K2 = — l/Aa7oi'2 along the boundary AlJq = {V2)~^^^ between the coherent domain and the AR2 
domain, whereas 7t is strongly reduced from K2 to Tci along the boundary AZOq = l/z^2 between the AR2 and ARl 
domains. 

Let us consider now the last case, corresponding to the domain V << 2 \ V2 {^^"^ for V1V2 < 0. In this case there 
is no spatial amplification in the coherent regime, but an absolute instability characterized by a spatial growth rate 
T^ab* ^ —1/ I ^2 1^/^. We thus find the same continuity as the one observed in the previous case V1V2 > along the 
boundary AcUo =| V2 \^^'^ between the spatial growth rate k"^*, corresponding to the coherent absolute instability, and 
the incoherent spatial am,plification growth rate K2; on the other hand k is strongly reduced from K2 to ki along the 
boundary between the AR2 and ARl domains, as in the case V1V2 > 0. We may thus conclude that in all cases the 
boundary AZUo = l/^'2, i-e. Au>o = 7o/j^2 in physical units, corresponds to a large reduction of the spatial amplification 
growth rate. 

Lastly, we may give for the sake of completeness the expression of the space growth rate k'^'"' associated with the 
incoherent absolute instability in the domain AR2. As stated before, for V1V2 < 0, the domain AR2 corresponds to a 

mixed situation where there is the coexistence of a spatial amplification characterized by K2 = 1 / AZUq | ^2 | , (for Vi > 
and V2 < 0) and of an absolute instability the space amplification growth rate of which is k"'"^ = 1/AcUo | V2 |^^^; one 
has thus I k"*"* |«| K2 \- 



4. Comparison of the Short Time and Long Time behaviors 

In Fig. 5 we reproduce the three boundaries found in the previous stability diagrams at which there is a reduction 

of the parametric growth due to the pump wave incoherence; the solid line corresponds to the boundary of the domain 
VR where the pump wave incoherence reduces the convective growth rate, i.e. where it modifies the short time 
behavior; the short- dashed line represents the boundary of the domain AR where the pump wave incoherence reduces 
the absolute growth rate, and the dotted-dashed line (which coincides with the solid line for F2 > 1) where it reduces 
the spatial amplification growth rate, i.e. the two domains where the long time behavior is modified. One immediately 
realizes that all the domains for which there is a growth reduction in the long time behavior are contained in the set 
for which there is a reduction of the initial growth, whereas the converse is not true. Thus the domain B in Fig. 5 is a 
domain where the spatial amplification is not affected by the pump wave incoherence although the convective growth 
is reduced in the initial stage; the same remark applies to the domain A for V1V2 > 0; for V1V2 < this latter domain 
corresponds to a regime where the absolute instability growth rate, although much smaller than the convective growth 
rate, is not modified by the pump wave incoherence whereas the initial convective growth rate is reduced. We may 
have therefore the surprising conclusion that the reduction of the parametric growth is more difBcult to achieve in 
the long time limit than in the initial stage. 
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5. Finite Length Thresholds for Absolute Modes: 

The thresholds set by the physical damping of the modes have been given for convective modes by (VI. 5) and for 
absolute modes by (VL9) in the incoherent limit or by (VII lb) in the coherent limit. In addition to satisfying these 
thresholds, the solutions must satisfy a threshold length criterion. 

In the case of absolute instability, the length criterion in the coherent limit is known to be 

tan ( ^,\7J°^" ) L = - ( ^ - 1 ) {VI. 29a) 



which in the absence of damping, 'jabs sets 



2 V'l 1^1 ^2 I 70. {VI. 29b) 



A similar condition is obtained by solving the set (V.32) with i/j = and with the boundary condition ni{l) = 
n2{L) = for ViV2,0 and Vi < 0. This RPA threshold condition is given by the unusual expression, 

27^ V Awi J V ^2 Aw2 / 

which, not unexpectedly, does not reduce to the coherent expression as Auij 0. In the limit of purely temporal 
incoherence, (VI.30) reduces to 



rRPA _ 1 r2_Aw Vi 



whereas, for purely spatial incoherence, (VI.30) reduces to 

Ai^L^ln I ^ I 
(I i I - I t I) 

The main point is that the ensemble averaged absolute growth rate for the infinite system given by VI. 8a) in the RPA 
theory can be obtained only if the plasma size L sufficiently exceeds a length given by the maximum of (VI. 29) and 
(VI.30). 



VII. DISTRIBUTION OF MODE AMPLITUDES AND INTENSITIES 

The statistical analysis has obtained dispersion relations from which growth rates, thresholds, and amplification 
rates in different domains of incoherence were obtained for the ensemble averaged behavior of parametric instabilities 
with incoherent pump waves. In this analysis, some questions arose such as the validity of the RPA dispersion relation 
in the intermediate domain that can be answered in principle at least by a direct numerical integration of the coupled 
mode equations (11-10). Not only do we test the assumptions and approximations of the statistical analysis but, in Sec. 
VIIA, we illustrate the meaning of the averaging procedure. The fundamental assumption that the RPA equations 
apply in the whole incoherent domain is examined in Section (VII. B); finally the relation of the statistical growth 
rates and thresholds to the growth rates and thresholds of a particular model of induced laser beam incoherence is 
examined in Section (VII. C). 



A. Ensemble Average Growth Rates and the Distribution of Intensities 

As has been noted at several points in this article, the incoherent absolute growth rates for the average intensity can 
be twice the rate expected on the basis of the growth rate for the average amplitude if certain inequalities are satisfied. 
For the purely temporal problem (where no distinction between absolute and convective is made), this factor of two 
arises if the damping is weak, specifically if Max {vj) < Jq/Aojq. However if this inequality is reversed, the average 
amplitude and average intensity growth rates are equal (we remind the reader that the definition of growth rates in 
Eq. III.9 removes the expected factor of two from the intensity growth rate). In the space-time problem considered 
in this article, one again notices a factor of two if | | = | V2 | in Eq. (VI. 8a) for purely temporal incoherence. 
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Aoji = /S.LO2 = AcjQ- If one realizes that these random processes lead to a distribution of mode amplitudes f(a) after 
a time T, then it is clear that the first moment (a) = J daf{a) and the second moment (| \) = j f{a) will be 
related only in special cases by (| |) =| (a) p. If the distribution remains narrow then these two moments are 
simply related. However we will discover that, in some cases, that is just when there is this factor of two at issue, 
the distribution is broad. The second moment then depends on finding the distribution of rare events where the 
wave amplitude grows to larger amplitude than the mean amplitude. We anticipate the result that the width of 
the distribution increases in time as does the mean so that the task of computing the ensemble average (| a? \) is 
increasingly more difficult as the sampling time is increased. 



1. Purely Temporal Evolution 



These remarks are now illustrated by considering the coupled mode equations for the purely temporal problem 
which have the simple form given by 

ai{t)=iaS{t)al{t), 



(^^ + a2{t) = 'yoS{t)al{t). {VII.l) 

In previous work, an exact solution for the ensemble average mode amplitude and intensity was found when S{t) is 
a Kubo- Anderson process (Section V.l). When Acjq > 7q Acjq > 1^1,1^2 the growth rates are found to be 7(a) and 
7((j2) = 2jq Aojq. Here, we are interested in the fluctuations of the amplitude and intensity about the means (| a |) 
and (I a p). Some useful insight can be gained into this question of fluctuations and its influence on the behaviors of 
I a I and | a p by considering the simplest possible model, that of two coupled undamped modes ai and 02 described 
by Eqs. (VII.l) where S{t) = exp{i(j){t)) and (j){t) is a random phase with bandwidth Auq- In this simple case we will 
be able to calculate the evolution of the distribution function f(a) of amplitudes. 
Because g =| ai p — | 02 P is a constant of the motion, we can conveniently write: 

ai = q^/^coshee''^' 



a2 = g^/^sinh^e*'^^ {VII.2) 
from which the equations of motion (VII.l) imply 

d9 

— = 70COSV' (^77.3) 



^ = -270 coth26' sinV' + ^ (VI I A) 
at at 

where 

^ = ^2 - (Ai + (/> {VII. 

If we define a distribution function f{9,tp,t), it evolves according to the Fokker-Planck equation: 

df d d 52 f 

+ ^7ocosV'/ - — 27ocoth2^sinV'/ = Awq^ {VI 1. 6) 

In the large bandwidth limit Awq >> 7o, d/dt., we can use a Chapman-Enskog like procedure, expanding / = 
/o + /i + ••• and to zero order obtain d'^fo/dip'^ = 0, so that /o = /o(^, i) independent of ip. To first order we obtain 

^"""'^ ~ 2coth20/o) (^77.7) 
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and to second order 

f - ^^r'''^^e {^0 - - |;27ocoth2.sin^A = (^77.8) 

Averaging over tp, we obtain the consistency equation, 

dfo _ 1 7o' d f d 
dt 2 Awo de 



coth26i /o (y77.9) 



We are interested in the evolution of f from an initial condition f{t = 0) = 6{9 — 9o)- If we assume 9 » 1, since 
we expect 9 to increase we can approximate coth2^ = 1 in the above equation. It is then readily solved obtaining 

/ = (27r7o't/Aa;o) exp (- - - 7oVAwo)' / (27o't/Aa;o)) ; {VI 1. 10) 

that is, 9 is Gaussian distributed with mean + lot/ ^<^Q and variance ^^t/ IS.u)q. This distribution is shown in Fig. 
7.1 for AcjQ = IO70 and 70^ = 50 together with some numerically generated distributions. Clearly then (a^^^) = 
^{0) = 7g/Aci;o and (a) ~ g^/^exp (^o + 7oV^'^o)- On the other hand, if we calculate (| a |") ~ g"/^(exp {n9)), we 
obtain 

(I a I") ~| a 1^ exp (n (n + 2) 7^t/2Awo) • (1^77.11) 

For the special case n = 2 we observe that the growth rate for (| a p) is 47g/Aa;o, containing the factor of two noted 
earlier. In particular wc see explicitly that the mean growth rate is definition dependent because the distribution of 
growth rates is so broad. To compare with these analytic results, Eqs. (VII. 1) were numerically integrated for a large 
number of different independent realizations of a Kubo- Anderson process such that l(a)T >> 1 and Awo = IO70. 
In Fig. 6, the distribution (broken line) of the intensity | a p is shown for T70 = 50 and 2 x 10^ realizations. The 
"measured" value of 7^^^ is 0.095 and 7(a2) is 0.16, and obviously the distribution at | a p= exp (47Qr/Aa;o) = 
exp (20) = 10*'^ is not resolved. Moreover, the few events with | a p> 10^'^ that have grown at the rate greater than 
0.165 70 have determined 7(a2)- In general, the accurate numerical computation of random multiplicative processes 
in which rare events determine the quantity of interest is nontrivial[77]. On the other hand, the distribution of events 
that have grown at the rate 7(a) is well resolved. When there is sufficient damping, i.e. Max (1^1,1^2) > IqI^^o the 
exact theoretical results show the ensemble average growth rate 7(a2> — 7(a)- The simulation results for V2 = 7o also 
shown in Fig. 6 (dashed-dotted line), demonstrate that the damping has dramatically narrowed the distribution of 
I a p so that 7(^2^ = 0.092 and f{o?) is zero for | |> 10^. 

Note that this behavior is not specific to the Kubo- Anderson process. Numerical simulations with a randomly 
diffusing phase or a K-A process with equally spaced elapsed times obtained similar distributions of | |. 



2. Space-time Evolution of Absolute Instabilities 

Here, the incoherent absolute growth rate 7(|a2|) for temporal incoherence is twice the "nominal" incoherent growth 
rate, 7i„c = lol if I ^1 1 = 1 ^2 |- Given the results just presented for the purely temporal problem, this factor of 
two suggests that we might expect a broad distribution to result from considering different realizations of a temporal 
random process with bandwidth Awq- In this case, we remind the reader that no absolute instability is allowed by the 
average amplitude equations since the average amplitude equations decouple in the RPA or Bourret approximation. 
It is also interesting to observe that, unlike the purely temporal case, damping on either mode does not affect the 
factor that multiplies to obtain 7(a2)- This factor, q~l + 'l^/VxV^I (| Vi | + | V2 |), depends on the ratio of the 
group velocities and approaches two when the group velocity ratio approaches one. Thus, in analogy with the purely 
temporal problem, we expect the distribution of intensities of growing modes to broaden as the group velocity ratio 
approaches one. 

Before we proceed to test this hypothesis, we must be certain that both the incoherent absolute growth rate 7(0^) 
(VI. 8a) is less than the coherent rate (Vl.IIa) and that the plasma length exceed sufficiently the threshold lengths 
(VI. 31) and (VI. 29a). In addition the threshold set by the physical damping and bandwidth must be exceeded. These 
conditions are met by choosing the parameters: Vj = 0, j = 1,2, Aujq = IO70, and L = lOLg. When the group 
velocity ratios V1/V2 = —16 and -100, the predicted incoherent absolute growth rates 7(|o2|) are 0.147 and 0.12 70, 
respectively. The predicted coherent absolute rates are 0.47 and 0.2 70 respectively. 

The distribution of | | for At^o = IO70 is shown in Fig. 7.2 for each velocity ratio. Each event in the distribution 
represents the ratio (T) =| ai {Xo,T) p / | ai {Xo,to) ^ where jo{T — to) = 45 and Xq is the initial position of a 
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source at t = in a2. The time to is chosen nonzero to avoid the influence of transients and was 20^^ in all these cases. 
This "initial" time is longer than the time for the convcctivc pulse to move across the system, t(, = 2i/ (| l^i | + | V2 |) 
but not always longer than the time it takes the slow wave to transit the length L. However since the growth is 
measured at the source point, the fact that in the entire plasma, the waves are not growing at the asymptotic rate 
should not matter. 

In these simulations, the numerically calculated ensemble average growth rate 7**™ was smaller than the RPA 

prediction in all cases as expected. When the velocity ratio V1/V2 = —16 and -100 the computed j^^T^ was 60% and 
75% respectively of the RPA prediction. However, we were most interested n the width of the distribution. When 
I V1/V2 I was 16 and 100, the full width at half maximum of the distribution was 7% and 2% of the peak respectively, 
and thus displayed the narrowing expected. When V1/V2 = —1, the computed distribution was not significantly 
broader than the distribution for V1/V2 = —16. Thus the difference in growth rates between V1/V2 = —1 and -16 may 
arise from differences in the unresolved tail of the distribution. 

We conclude from this study that the RPA growth rates for the ensemble average intensity are correct but in a 
given simulation, or perhaps a particular experiment, the measured growth rate may well be a factor of two smaller. 
However when the group velocity is large, as is typically the case of physical interest, the RPA predictions are more 
representative of the typical behavior. 

B. The RPA Conjecture 

The question here is whether the RPA equations are applicable in the whole domain defined by the inequalities 
Auii > Max(7o, i^i) and Au>2 > Max (70, 1^2) or whether there is an intermediate domain when one of the cross 
inequalities Awi > 1^2 or AuJ2 > vi is not satisfied. To examined this case we look at the spatial Kubo-Anderson 
(K-A) process analyzed previously^ with zero temporal incoherence and spatial correlation function with the property 

{S {x) S* {x + Ax))^ = exp {-AKAx) . {VII.12) 

If we assume | Vi \»\ V2 \ and V1V2 < 0, then the intermediate domain is described by the chain of inequalities, 

AKVi > vi >\ AKV2 \> z^2,7o- (F//.13) 

A physical example that might produce this set of inequalities is two plasmon decay in the presence of stationary ion 
acoustic turbulence where the shorter wavelength Langmuir wave is more strongly damped and has a faster group 
velocity. The choices in our numerical solution were: 

AK I V1V2 /70 =4, z/i = 270, 1^2 = 0, V1/V2 = -16, [VII.U) 

for which the RPA equation absolute instability growth rate 7(a2> = O.I270 and the absolute instability coherent 
growth rate (VI.ll), 7°^^ = 0.8670. The simulation region L = 8Lg, [Lg =\ V1V2 |^/^ /70) was larger than the 
threshold length given by the maximum of Eq. VI. 32 and VI. 29. The simulation was run until 7ot = 50 for 1000 
independent realizations of this spatial K-A process. An initial value was given to the undamped wave at t = 
but the intensity at jqT = 50 was measured relative to the intensity at joto = 20 because by that time the modes 
were observed to grow at this time asymptotic rate. The ensemble averaged rate 7{a2) was found to be independent 
of whether its calculation was based on the intensity at the source point | ai {Xo,T) ^ or the total mode "energy" 
J dx \ ai {x,T) p. The results support our conjecture that the RPA equations apply in the whole incoherent domain, 
because the measured intensity growth rate, 7^*"2" = O.II70 agreed with the RPA rate J{a^) ~ O.I270. However, the 
distribution of intensity at 70T = 50 (relative to 7oto = 20) displayed in Fig 8 showed that, in fully a third of the 
cases, no growth occurred and that the remaining distribution consisted of a slowly decreasing tail. The maximum 
intensity in the distribution was achieved by a mode that grew at close to the coherent rate. In addition, the faster 
growth rates were achieved in cases in which there was an abnormally long distance between phase changes; in the 
maximum growth case, this distance was 2.18 Lg. As we remarked at the outset of this section and in the introduction 
(Sec. I), the growth of parametric instabilities with this model of spatial incoherence was studied by Williams et al. 
using more sophisticated statistical methods (quite different in nature to those used in this article) with the result 
that the fastest growing mode grows at the rate 7/, given by Eq. 1-9. This rate 7/ depends on the size of the system, 
which is understood on the basis that as the system gets larger there is an increasing possibility that there will be 
a large region with no phase change. In our simulation this size dependent factor InAKL = 3.47 turns out to be 
approximately equal to the numerical factor of 4 in the RPA growth rate (Eq. VI. 8a); all other factors are identical 
and AK Ako . Unfortunately it is not feasible to simulate a system that doubles the size of this logarithmic factor 



34 



with a constant value of Aif. Thus these simulations do not distinguish between the RPA theory and the results 
of Williams et al. However, the simulations do support the conjecture that the incoherent results apply in whole 
incoherent domain. 

C. Modeling of Induced Spatial Incoherence 

One motivation for this work is the induced spatial incoherence technique (ISI) for creating smoothed laser beam 

intensity distributions in the focal plane. In this technique, statistically independent beamlcts overlap and thereby 
produce both intensity and phase variations. Only pumps with phase variation alone have been considered up to this 
point in our numerical examples. Here we wish to show that the same statistical methods are valid for intensity and 
phase varying pumps if the conditions for the incoherent results to apply are met. 

A model that exhibits both intensity and phase variation in one dimension is given by, 



when 5Kj is uniformly distributed on and (j)j{t) is a random variable of time. Each (t)i varies independently 

with the same mean coherence time, (Awq)"^- The spatial coherence length in this model, Xc = ir/Km, is related to 
the focal length of the focussing optic. A shorter focal length results in a shorter coherence length. 

At any given point and instant of time, the intensity | 5* (.t, t) being the absolute square of the sum of N complex 
numbers of length N~^/^, is Gaussian distributed with unit mean and variance. Because the sum varies with position 
over a length scale x^, the instantaneous spatial intensity pattern is quite spikey with peak to average variations of 
more than two quite common. However because the sum at a given spatial position also varies in time with coherence 
time (Awo) ^, the time averaged intensity is a much smoother function of space. Clearly if (t) is constant in 
time, then the spikey intensity pattern of Eq. VH.IS is stationary which can lead in the coherent domain to growth 
enhanced relative to a uniform intensity. Even in the incoherent domain without temporal bandwidth, we saw in the 
previous section that there is a finite probability, especially in large systems, for the phase to remain constant for a 
large distance. With phase and intensity variations, large fluctuations in growth rates are expected in the absence of 
temporal bandwidth. In the complementary limit where Aui > A0J2 > Max (70, 1'l, ^2), the results of our statistical 
analysis should apply. For simplicity, we choose to compute the convective growth for both temporal and spatial 
incoherence in the incoherent limit and for large group velocity ratios, | V1/V2 |>> 1. With a single term in the 
sum in Eq. VII. 15 and Awi = Al02 = Acoq, the convective growth rate was in agreement with Eq. (VI. 7). With no 
temporal bandwidth, Awq = 0, the pump intensity has striking amplitude variation with peak to average variations 
of more than three to one. Nonetheless if | AKV2 |> 7oj the convective growth rate, is reduced according to Eq. 
(VI. 7) even if Awq = 0. As the spatial incoherence weakens and the boundary AKV2 = 7o is approached for Awq = 
the measured growth rate in the simulations diverges from Eq. (VI. 7) until, for 7q >| AKV2 |, growth rates in excess 
of 7q can occur with large fluctuations from case to case. Laser bandwidth, even if Aluq < 7q, reduces both the 
fluctuations and the growth rate to less than 7q. Thus, some bandwidth is needed in this model to provide a smooth 
transition from the coherent to incoherent domain. Other models such as the one used in Section VII. B have the 
property that | S {x,t) |= 1 so that the statement that pump incoherence can only reduce the growth rate is indeed 
true. Simulations in that case do show a smooth transition from the coherent to incoherent growth rate 

VIII. CONCLUSIONS AND APPLICATIONS 

A comprehensive treatment of the effects of pump wave temporal and spatial incoherence has been presented for a 
homogeneous plasma. Eschewing specific models of incoherence that allow exact solutions, we have derived equations 
for ensemble averaged mode intensity and amplitudes in the incoherent limit. At each step in the process of deriving 
these approximate statistical equations, the conditions that must be satisfied are clearly stated. Thus, the meaning of 
incoherence in the context of paranic;tric instability theory is carefully defined. The set of inequalities that comprise 
this definition form a major contribution of this work. Of course, the primary contribution of this work is the set 
of thresholds, growth rates, and amplification rates for both coherent and incoherent, convectively and absolutely 
unstable parametric interactions. 

However, the majority of our readers are interested in the more practical results concerning thresholds, growth 
rates, and amplification coeflacients as they pertain to laser plasma interactions. In a subsequent work,^ we will apply 
the results presented here to particular instabilities, and include both backscatter and sidescatter geometries. Here, 
we briefiy consider the examples of stimulated Raman and Brillouin backscatter for parameters of interest to laser 




(^77.5) 
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fusion. Wc choose the electron temperature, Tg = 3keV, the charge state Z = 40, the atomic number A = 80, the 
laser wavelength, Ag — 0.35/xm, the laser intensity I = 10^^ W/cm^ and the electron density, Ng = 0.2N(. for SRS and 
0.25 Nc for SBS. Here Nc is the critical electron density defined by A-kNcB^ /me = u)q. For almost all calculations it is 
the ratio Z/A and not Z or A by themselves that is important. 

The nominal SBS growth rate is 7q''* = 2.5 x 10~^u;o whereas the SRS coherent convective threshold, Eq. (I-l) 
is = 3.2 X 10~'*wo- Landau damping was neglected because for Te = 3keV, an electron density can always be 

chosen to make it negligible. Thus SRS is an order of magnitude above threshold. The laser bandwidth sufficient to 
achieve convective stability can be found by using Eq. (VI.5). We find AtOc = 4.5 x 10~^a;o- On the other hand, 
the coherent threshold for absolute instability given by Eq. (VI. 9a), 7oahs = 1.510^^ujq, is exceeded by less than a 
factor of two. However, as we have emphasized, in Sec. VI, relatively large laser bandwidth is required to reduce 
the absolute growth rate. Satisfying the validity conditions for using the incoherent formula, Eq. (VI. 10), imposes 
the condition, Aluq/luo = 10~^, which is, coincidentally, the same condition numerically that is required for absolute 
stability as determined by Eq. (VI. 9a). 

The SBS thresholds are dependent on the ion acoustic damping rate and thus the ion temperature, Tj. However, 
even if Tj = Tg the ion Landau damping is negligible for Z > 10. Thus, the ion acoustic damping rate is determined 
by the electrons and is given by 7oc/wo = 2.3 x 10~^. Hence, the convective coherent threshold joconv/i^o ~ 1-5 x 10~^ 
is easily exceeded by the nominal SBS growth rate 70 = 7 x IO^'^wq- In addition, the SBS absolute instability 
coherent threshold 7oo6s = 3.4 x 10~''a;o is exceeded. Here, the bandwidth required to achieve convective stability is 
AiVc/oJo = 2.2 X 10~^. Moreover, the bandwidth required for absolute stability is by accidental choice of parameters, 
the same as for convective stability because although i>i << V2, nonetheless {V1/V2) 1^2 > and the incoherent 
absolute threshold Eq. (VI. 9a) reduces to the incoherent convective threshold Eq. (VI. 6). At this point, we should 
warn the reader that the results reported here assumed the validity of the coupled mode equations which requires the 
mode frequencies Wj > Aujq for j = 0,1,2. However, we have in the case of SBS a critical bandwidth for stability that 
exceeds the ion acoustic frequency = 2 x IO^'^luq. At the very least, one might expect additional stability when 
AojQ > oJa- For a partial answer to the effects that Aujq > uja has, we considered SBS again for the purely temporal 
problem without using the assumption that u)a » Auq. We start with the equations. 



^ + a2 = -27oWaat, {Vlll.lb) 

where 5uj = ojq — oJi, is the difference between the pump frequency and the scattered light frequency. Here ai,2 are 
the mode amplitudes for the scattered light and ion acoustic waves respectively. Two different dispersion relations 
are obtained, as expected, by obtaining equations for the evolution of average amplitude (aj). The ion acoustic mode 
(02) dispersion relation is not affected by magnitude of Awo/^a so that the growth rate, 

^<-) = -'^^+(^;rTk)' ^""'''-'^ 

is a maximum when 6u) = oJa- On the other hand, the dispersion relation for (ai) is modified so that the growth rate 
when Awo >> Wa, f2, 

3V3wa7o /T/^■r^-Q^ 
7(ai> = + — A^' iVIII.3) 







shows a stronger reduction with bandwidth than 7(^2) . Here the maximum growth rate for (ai) occurs for a bandwidth 
dependent real frequency — = Awq/VS- Thus, the threshold for mode (ai) is increased when Awq > by the 
ratio Auio/iVa- The actual threshold for the coupled set VIII. 1 is, 

7o'con. = Min (^U2 {i^i + Aojo) , lyi {V2 + Awo) Max (\, ^^^i^±^^ ^ , {yillA) 

This equation (VIII. 4) is a generalization of Eq. (VI. 6). In our example, the threshold is not increased when Ao^o > ^a- 
In our example, chosen for its relevance to laser fusion, the critical bandwidth to reach convective or absolute 

stability has been found to be 2-5% of the laser frequency. Such a bandwidth is beyond the currently achievable with 
any ICF laser system. However, novel techniques are being pursued that are capable of achieving these bandwidths[6]. 
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In addition, there is the question of whether the average laser intensity or the larger hotspot intensity should be used 
in evaluating the thresholds. The use of the larger intensity would make the bandwidth requirements infeasible. 
However, if bandwidth is combined with a beam smoothing scheme, use of the mean intensity may be correct. Recent 
calculations[66,67] have shown that, for at least one such scheme, ISI, the use of the average intensity is justified. 

In addition to temporal bandwidth, there is also the effect of spatial incoherence produced by beam smoothing 
techniques that lead to a spread of wavevectors, Ako- In the examples of backscattcr considered previously in this 
section for temporal bandwidth, spatial incoherence will have little effect because Afco**/Afco* ~ /*~^ where the f 
number of the focussing optic, f#, is typically large in ICF applications. Furthermore, the wavenumber spread will 
not increase the convective threshold because the effective bandwidth, | Ako ■ Igmin |, is in general small compared 
to 1^1,2^2) and 70- However absolute instabilities with significant wavevectors perpendicular to the laser propagation 
axis such as two plasmon decay may be influenced significantly by such a spread since there Afc*V^*i„ax will be the 
effective bandwidth. 

The work reported here is only a step in the application of a statistical theory to particular instabilities. Further 
progress may require modifications or extensions to our work in order to remove some of the limitations of our 
approximations. We look forward to new developments in the theory and even more to the stimulus of experimental 
data. 
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FIG. 1 Geometry of the scattering: Kq, Ki, K2 are the mean wavenumbors of the pump wave and of the decay waves; the 
pump wave vector Kq varies in magnitude an amount A | feo | and lies in a cone of angle A^o about the direction of the mean, 
Ko ■ Also shown are the angles a2 and ai which the mean wavevectors Ki and K2 make with the mean pump wavevector Kq 
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FIG. 2 Schematic stability diagram for the convective growth rate. The sohd hne corresponds to the boundary of the domain 
VR where the pump wave incoherence reduces the convective growth rate, and the dotted hnc is the boundary of the domain 
VS where it completely stabilizes the parametric coupling. The rightmost vcrticlc line divides the coherently stable (CS) from 
the coherently unstable (CU) domains. This diagram assumes vi > vi corresponding to the ordering assumed in the main text. 
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FIG. 3 Schematic stability diagram for the absolute growth rate. The short-dashed line, between I V on the vertical axis 

and 2 I w 1^/^ on the horizontal axis, represents the boundary of the domain AR where the pump wave incoherence reduces the 
absolute growth rate. Here v =\v2 \- The long-dashed line is the boundary of the domain AS where the pump wave incoherence 
completely stabilizes the absolute growth. To the right of the rightmost vertical line is the domain of coherent stability. This 
diagram assumes 172 >V\. 
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FIG. 4 Schematic stability diagram for the spatial amplification growth rate. The dottcd-dashed lines represent the boundaries 
of the domains ARl and AR2 where the pump wave incoherence reduces the spatial amplification rate; the dotted line of crosses 
is the boundary of the domain AS (identical to the domain VS of Fig. 2) where it completely stabilizes spatial amplification. 
The short-dashed line (see Fig. 3) represents the threshold for the existence of absolute instabilities in the case V1V2 < 0. The 
rightmost vertical line separates the coherently stable (CS) and unstable (CU) domains. This diagram assumes V2 >V\. 
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FIG. 5 Comparison of the short time and long time behavior. The domains A and B correspond to regimes for which the pump 
wave incoherence reduces the initial parametric growth (convective growth) whereas the long time behavior remains unaffected. 
The meaning of the solid and dotted, the short-dashed, and the dotted-dashed lines are given in Figs. 2, 3, and 4 respectively. 
The diagram assumes F2 >V\. 
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FIG. 6 Distribution of the intensity | a\ ^ for the purely temporal problem. The solid curve is the analytic solution, and the 
broken curve is the numerically generated solution for undamped waves. The dashed-dotted curve is the numerical solution 
when V2 = 7o- The bandwidth is Awo = IO70 and the abscissa is the base 10 logarithm of the intensity | ai |^ after an interval 
7oAt = 50. 
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FIG. 7 Numerically generated distribution of the intensity | ai p for the space-time problem with purely temporal bandwidth 
Awo = IO70 after an interval 70 At = 45. The solid curve has | V1/V2 |= 16 and the dashed curve has | V1/V2 \= 100. The 
abscissa is the base 10 logarithm of the intensity | oi f. 




FIG. 8 Numerically generated distribution of the intensity | ai p for a spatially incoherent pump with Afco | yj\ V1V2 I/70 = 
4, vi = 270, 1^2 = 0, and V1/V2 = —16. The abscissa is the base 10 logarithm of the intensity | ai \^ after an interval yoAt = 30. 



